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SPREADING ESTIMATES FOR QUANTUM WALKS ON THE 
INTEGER LATTIGE VIA POWER-LAW BOUNDS ON 
TRANSFER MATRICES 

DAVID DAMANIK, JAKE FILLMAN, AND DARREN C. ONG 


Abstract. We discuss spreading estimates for dynamical systems given by the 
iteration of an extended CMV matrix. Using a connection due to Cantero- 
Griinbaum-Moral-Velazquez, this enables us to study spreading rates for 
quantum walks in one spatial dimension. We prove several general results 
which establish quantitative upper and lower bounds on the spreading of a 
quantum walk in terms of estimates on a pair of associated matrix cocycles. 
To demonstrate the power and utility of these methods, we apply them to 
several concrete cases of interest. In the case where the coins are distributed 
according to an element of the Fibonacci subshift, we are able to rather com¬ 
pletely describe the dynamics in a particular asymptotic regime. As a pleasant 
consequence, this supplies the first concrete example of a quantum walk with 
anomalous transport, to the best of our knowledge. We also prove ballistic 
transport for a quantum walk whose coins are periodically distributed. 
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1. Introduction 

In recent years, quantum walks have been studied extensively; see [2 El 0 0 H 
iniiisiiiiiiiiisniiniiMiEiissiiiiiin] for some papers on this subject that have 
appeared in the past five years. These are quantum analogues of classical random 
walks. For simplicity, let us focus on the important special case of walks on the 
integer lattice Z. 

A classical (stationary, or time-homogeneous) random walk on the integer lattice 
Z is given by transition probabilities {pn,m}n,mGZ, where p„^rn S [0,1] denotes the 
probability of moving from site m to site n. Denote the bi-infinite matrix with these 
entries by P. We therefore require that each column sum of P is equal to one, that 
= 1 for every m G Z. The state of the system may be specified by a 
column vector v = {vn)n& with entries Vn G [0,1], which is normalized in the sense 
that X^Tiez iin = 1 so that Vn describes the probability of the system being at site 
n. The time evolution of this state is then given by v{k) = P^v for A: > 0 or, if P 
is invertible, for A: G Z. 

An important special case is obtained by only considering nearest-neighbor tran¬ 
sitions with symmetric transition probabilities. That is, we have Pn,m = 0 for 
|n — m| > 1 and Pn,m = Pm,n- In this case the matrix P is tridiagonal and symmet¬ 
ric and hence has the structure of a Jacobi matrix. To einphasize this fact we will 
write in this case J instead of P for the transition matrixo Thus, in this case, the 
evolution of an initial state v is given by v{k) = J^v. In particular, the probability 
of being at site n at time k is given by 

( 1 . 1 ) 

where dn is the vector with value 1 at position n and value 0 otherwise, and 
{v^,v^) = 'Yhn&'^h'^n- Here is where our symmetry assumption turns out to be 
critical. Since J is obviously a bounded self-adjoint operator on ^^(Z), an expres¬ 
sion like dm can be rewritten with the help of the spectral theorem. Indeed, it 
takes the form dp,n,j,v{x), where Pn,j,v is a complex measure supported by the 
spectrum of J. Therefore, spectral theory enters the game and the long-time evo¬ 
lution of the classical random walk can be studied by means of a spectral analysis 
of the transition operator J. 

To motivate the discussion of the quantum analogue, we can interpret the tran¬ 
sitions as follows. Assuming further that the diagonal elements are zero as well, we 
have Pn-i,n + Pn+i,n = 1- Thus the transition from n to one of its neighbors nil 
corresponds to tossing a coin with probabilities pn-i,n and 1 —pn-i,n = Pn+i,n and 
transitioning according to the outcome of this coin toss. A quantum walk on the 
integer lattice is described by assigning a spin to every site n G Z, which is given by 
a vector, and determining the transitions from n to one of its nearest neighbors 
nil with the help of a “quantum coin,” which is a unitary 2x2 matrix. The 
details are described in Subsection 12.41 below. The upshot, following an important 
observation of Cantero, Griinbaum, Moral, and Velazquez in 0, is the following. 
The time-evolution may now be described by an extended CMV matrix £ in place 
of the Jacobi matrix J above. CMV matrices are the natural unitary analogues of 

^For a Jacobi matrix, one usually requires in addition that the off-diagonal terms are strictly 
positive. This would be in fact a natural additional assumption since the case an := Pn,n—l = 0 
will correspond to two non-interacting walks on Z n (—oo,n — 1] and Zn [n, oo), which can be 
studied separately. 
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Jacobi matrices, and they play a canonical role within the class of unitary opera¬ 
tors analogous to the canonical role played by Jacobi matrices within the class of 
(bounded) self-adjoint operators; see [331 HI] and references therein. Thus, we are 
now concerned with the study of 

(1.2) (<5„,fS), 

where £ is an extended CMV matrix, and in particular a unitary operator on 
Again, this allows one to employ spectral theoretical methods to analyze 
the behavior of the corresponding quantum walk. In particular, one may relate 
continuity properties of spectral measures with respect to Hausdorff measures on 
the unit circle 9D to spreading properties of the evolution (II.2|) and hence of the 
quantum walk; this approach was developed in m and applied to concrete quantum 

walks in m- 

The approach just mentioned works for arbitrary unitary operators. For ex¬ 
tended CMV matrices, one can sometimes do better and obtain improved spread¬ 
ing estimates by using a different method. While the approach above is based on 
rewriting dOJ) by means of integrals over the unit circle with respect to spectral 
measures, the alternative method we propose and develop in this paper rewrites 
(EH) by means of integrals over the unit circle with respect to normalized Lebesgue 
measure (i.e., the measure generated by normalized arc length). The fundamental 
formula is given in Lemma l3.II below and was already observed in Eia It connects 
suitable averages of the quantities (11.21) to integrals involving matrix elements of 
the resolvent of £. This in turn connects time-averaged spreading to growth prop¬ 
erties of transfer matrices, as the latter matrices are related to the properties of the 
resolvent probed by the integral in the fundamental formula. This connection is 
exceedingly useful in many cases of interest as there is a plethora of tools one can 
use to study the growth of transfer matrix norms. These tools are CMV analogues 
of tools originally developed in the Jacobi setting. In fact, some of these CMV 
versions are developed in the present paper. 

This paper consists of two parts. The first part develops the general theory 
which enables one to deduce spreading estimates for the discrete time dynamics 
generated by an extended CMV matrix from bounds on transfer matrix norms 
for spectral parameters off the unit circle, which immediately gives general results 
when specializing to the case of quantum walks on the line (precise descriptions 
of the results may be found in Section jj). These general results are analogues for 
extended CMV matrices of results known for discrete Schrbdinger operators. The 
latter operators generate a continuous time evolution via the Schrodinger equation. 
The connection between discrete Schrodinger operators, or more generally Jacobi 
matrices, and CMV matrices is very fruitful and has been explored extensively in 
recent years. Most of the time results are proved first in the self-adjoint context and 
are then carried over to the unitary context. There are some notable exceptions, 
such as Rakhmanov’s theorem. Carrying over results from the self-adjoint case 
to the unitary case, or more concretely from the Jacobi to the CMV setting, is 
sometimes straightforward, sometimes prohibitively difficult, and sometimes doable 
but far from straightforward. What we do in this first part of the paper falls in 


^The formula was proved in m in the hope that it would eventually become useful once one 
is able to work out the unitary analogues of a series of papers by Damanik and Tcheremchantsev 
in the self-adjoint case. The present paper realizes this vision. 



4 


DAVID DAMANIK, JAKE FILLMAN, AND DARREN C. ONG 


the middle ground of doable but challenging. We obtain complete CMV analogues 
of the known Schrodinger results, but proving them required us to overcome quite 
a few obstacles. This is one of the cases where the ability to carry over the result 
should not be taken for granted but rather appreciated. 

The second part presents various applications of the general theory. These ap¬ 
plications range from simple but interesting observations to results that require a 
significant amount of work. Namely, in Sections |6] and [H we will discuss lower 
bounds for the spreading of quantum walks associated with polymer models and 
the Thue-Morse substitution sequence. These results are relatively easy to estab¬ 
lish and rely on the presence of a small set of exceptional spectral parameters at 
which the required estimates for the transfer matrix norms hold and are straight¬ 
forward to verify. Nevertheless they are interesting because the spectral measures 
are either difficult to study (in the Thue-Morse case) or they are so singular that 
the approach from m gives no non-trivial consequence for the spreading rates (for 
typical realizations of random polymer models). These examples demonstrate the 
usefulness of the method developed in this paper, which does not rely on spectral 
continuity properties and yet may be able to produce pretty strong estimates. 

The applications that require a substantial amount of work consider the case 
of quantum walks with Fibonacci coins or with periodic coins. Quantum walks 
with coins generated by a Fibonacci substitution sequence, discussed in Section [71 
are interesting because they can exhibit anomalous transport in the sense that the 
transport exponents (to be defined in the next section) take fractional values. To 
show this, one needs to establish upper and lower bounds on the spreading behavior 
of the quantum walk. This is another important advantage of the method developed 
here over the spectral continuity approach. The latter is incapable of establishing 
upper bounds purely in terms of spectral continuity (the case of random polymer 
models neatly demonstrates this). The method put forward here, however, natu¬ 
rally can be used to produce two-sided estimates, and in the Fibonacci case we work 
out in detail how this may be implemented. This requires quite extensive analysis 
(which is of dynamical and combinatorial nature) but in the end it produces rather 
sharp estimates, especially in the regime studied here. Our results rigorously estab¬ 
lish transport behavior which was predicted numerically in the physics literature 

m- 

Finally, the case of quantum walks with periodic coins is clearly a fundamental 
one, for which one expects ballistic behavior to occur. That is, the walk spreads 
out in space with a well-defined strictly positive velocity. In Section [51 we work out 
a result of this kind for periodic extended CMV matrices and derive the statement 
for quantum walks with periodic coins as a corollary0 While this happens only 
in the very last section of the paper, we want to emphasize that the proof is in 
fact independent of the bulk of the paper. Readers familiar with extended CMV 
matrices and the CGMV connection between quantum walks on the integer lattice 
and extended CMV matrices may skip ahead to Section [Q] 


^During the review of this paper, it was kindly pointed out to us that the statement about 
ballistic transport for quantum walks on the line with periodic coins is a particular consequence 
of (T] Theorem 4]. While the two proofs are related, we have decided to still present ours, which 
is inspired by and [E], in this paper for the convenience of the reader since the authors of [I] 
work in a more general setting and their paper does not state some particular consequences of 
their Theorem 4, such as our Theorem 19.11 and its corollaries in Section O 
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We conclude the introduction with a brief remark on the half-line case. One¬ 
sided, or standard, CMV matrices may be used to study quantum walks on Z+, 
as shown also in [^. The general theory we develop in part one of the paper has 
a natural half-line counterpart. In fact, some of the proofs are slightly simpler in 
this case. On the other hand, the applications we discuss in the second part of 
the paper all belong to the family of dynamically dehned operators over invertible 
dynamics, and therefore are more naturally studied in the two-sided setting. One 
could study their half-line restrictions, but one would encounter possible discrete 
spectrum which changes the dynamical picture for initial states from the associated 
spectral subspace. To avoid cluttering the paper, we state and prove all theorems 
only in the two-sided case, and we describe how the theorems may be adapted to 
the one-sided setting. 
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Part I. General Results 


2. Precise Statements 


2.1. Unitary Dynamics. We shall adopt the notation and development found in 
mm- Consider a unitary operator U acting on the separable Hilbert space "H, 
equipped with an orthonormal basis {(/^nlnGZ- In this paper, we will primarily focus 
our attention on two closely related scenarios: 

(1) t/ is an extended CMV matrix acting on 'H = equipped with the 

standard orthonormal basis ipn = Sn, n € h. 

(2) U is the one-step update rule of a quantum walk on Z, and TL = £^(Z)(g)C^. 
If we denote the standard basis of by {e-|-,e4,}, then the elementary 
tensors of the form 

(2.1) ip2m = Sm® e-f, (P2m+1 = <5m ® 64 ,, m £ Z 

comprise an orthonormal basis for 'H. 

Given an intial state ip €% normalized by ||'!/'|| = 1, we are interested in the time 
evolution of the vector ip, that is, we want to study the evolution of ip{k) = U^ip 
as fc £ Z+ grows. To quantify this, we first put 


a^{n,k) = \{ifn,ip{k))\'^ , k G Z, n G Z, 


which can be thought of as the probability that ip is in the state ipn at time k. We 
shall also be interested in the time-averaged probabilities, given by 



( 2 . 2 ) 


6 


DAVID DAMANIK, JAKE FILLMAN, AND DARREN C. ONG 


Taking time averages of quantum dynamical quantities has a long tradition; com¬ 
pare, for example, the survey-type papers [8j [20l [36] and references therein. On a 
technical level, this has its roots in Wiener’s theorem. In particular in the case of 
singular continuous spectral measures, Wiener’s theorem allows one to prove rig¬ 
orous estimates for time-averaged quantities, whereas in general similar estimates 
may not hold for non-time-averaged quantities. Additionally, the process of averag¬ 
ing can transform dynamical quantities in particularly pleasant ways. For example, 
the exponential averages considered in this paper recast the a’s into averages of 
the Poisson kernel against spectral measures; this point of view is explored in more 
detail in [9|. 

To be completely proper, the normalizing factor in (12.2|) should really be 
(l-e-2/^) so that a is also a probability distribution on Z. Of course, 2/AT 
is correct to first order in AT, and we are primarily interested in quantities in the 
limit A' —>■ oo, so this distinction does not really matter. In particular, this does not 
affect the /3’s (which will be defined below). We will be interested in such averages 
throughout the paper, so, for a function / : Z>o —?► K., we introduce the notation 
(/) to denote the average of /. More precisely, we set 

k=0 

For example, in this notation, one could write 'a^{n,K) = {a^[n, ■)){K). 

Since C/ is a unitary operator, = I for every k. We may then think of 

as defining a probability distribution on Z. From this point of view, it makes 
sense to describe the spreading of these distributions in terms of their moments. 
More precisely, for p > 0, fc € Z, and K G Z+, define 

IXI^4k) = J^(lnr + l)a^(n,k), 

(lXi;)(K) = J2(lnr + l)a^(n,K). 

We are interested in the spreading of states which are initially well-localized in the 
sense that |X|^(0) < -l-oo for all p > 0. Consequently, we define 

S = S(?{) := I€ ?{ : lim |n|^|(ip„, V))! = 0 for all p > 0 

|n|—>-oo 

Obviously, S depends on the choice of basis, but we will not vary the orthonormal 
bases of £^(Z) and £^(Z) 0 in this paper, so we will not make this dependence 
explicit in the notation. 

Remark. It is helpful to observe that the probability of finding the wave packet 
outside of a ball of radius R can be used to find a simple lower bound on the 
moments, viz. 

(2.3) \X\l(k) >R-Y: a.,p{n, k), Vi? > 0, fc G Z, 

|n|>R 

(\Xf^){K)>RP Mn,K), 

\n\>R 


(2.4) 


Vi? > 0, a: G z+. 
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In the two situations described above, there is a universal ballistic bound on the 
associated unitary dynamics; compare [201 Theorem 2.22]. More precisely, if U is 
an extended CMV matrix or the update rule of a quantum walk, ij) £ S, and p > 0, 
there is a constant C = C'^,^ > 0 such that 


1^1^ (fc) < CkP for all fc e Z+. 

Consequently, we would like to compare the growth of the pth moment to polyno¬ 
mial growth of the form for a suitable exponent /3 € [0,1]. In light of this, the 
following transport exponentt0 are natural objects to consider 


^^(P) = limsup 

k—^oc 

Pi, (P) = 1™ inf 

Plip) = limsup 
K—^oo 

Pi,(p) = liminf 

^ K —Vnn 


log 


plog(A:) 

log ( 

Vi;(fc)) 

plog(/c) 

log 

((ai:;)(m) 


piog(A:) 

log( 

<*■>) 


A-s-oo p\og{K) 


By Jensen’s inequality, /3^ and (3^ are all non-decreasing functions of p [20l 
Lemma 2.7]. 


2.2. CMV Matrices and the Szego Cocycle. Given a sequence {anjnez of 
complex numbers where anSD = {zGC: \z\ <l}for every n G Z, the associated 
extended CMV matrix, £ = is a unitary operator on £^(Z) defined by the matrix 
representation 


( . 


(2.5) £ = 


a?oP-i 

—Q0Q_1 oTpo 

PiPo 





PoP-i 

— Poo_i — oToo 

—plQO 






02P1 

— 0201 

03P2 

P 3 p 2 




p2pl 

— P2Q1 

—75iEa2 

—P 30 L 2 






CMp 3 

—0403 

OtSpi 

P 5 P 4 




P 4 P 3 

— P 4 Q !3 

— 0504 —phOtA 

V 







= (1- 

for every n 

G Z. 

Since we 

heavily 

use the 

paper [3^ 


where Pn = (t — ' lor every n ^ iL. Since we neavily use tire paper m 

what follows, we adopt their convention for the location of the diagonal elements, 
that is, we have {5n,£5n) = —otn+iOin for n G Z (so that our is their —a7()- 
We refer to {onlnez as the sequence of Verblunsky coefficients of £. Note that all 
unspecified matrix entries are implicitly assumed to be zero. As we have said in 
the introduction, one may also study (half-line) CMV matrices, which are given by 
setting ai = —1 and restricting £ to £^(Z+). 


^Some authors consider Cesaro averages for the moments, instead of the exponential averages 
which we consider. However, it is not hard to see that either method of averaging yields the same 
values for /3^; compare [201 Lemma 2.19] 
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For a S D and z S 915, the corresponding Szegd matrix is defined by 

= P=V'1-|«P- 

The Szegd transfer matrices associated to Ea are then defined by 


S'(a„_i,z)-' 

■ * S{am:Z) 

n > m 

/ 


n = m 

S{an,z)~^ ■ 

S(^CXrn—l ■) z) 

n < m 


where n,m € Z. These matrices arise in the setting of half-line CMV matrices as 
propagation matrices of the orthonormal polynomials associated to the (5o spectral 
measure. That is to say, Tin, m; z) may be used to obtain the degree-n orthonormal 
polynomials from the degree-m orthonormal polynomials; see |431144) for a much 
more comprehensive discussion. The Szegd transfer matrices are still quite relevant 
to the spectral analysis of extended CMV matrices; for example, the spectrum of 
an extended CMV matrix may be characterized via exponential growth properties 
of the Szegd matrices P^ . 

Our first main result says that quantitative polynomial bounds on the Szegd ma¬ 
trices imply quantitative lower bounds on the spreading of associated wave packets 
defined by iterating £. The precise formulation follows. 


Theorem 2.1. Let £ = Eq, be an extended CMV matrix whose Verhlunsky coeffi- 
eients are bounded away from 9D in the sense that 

(2.6) M := Halloo < 1. 

Suppose further that there exist C > 0, 7 > 0, such that the following condition 
holds: for every R>1, there exists a nonempty Borel set An C 9D sueh that 

(2.7) sup max ||T(n, m; z)|| < C'i?'''. 

z^Ar |"|.|m|<R 

For each K G Z_|_, let Rk = and let Bk denote the 1 /K-neighborhood of 

the set Ajir in 915, i.e., 

Bk = G 9D : |z — z'I < ^ for some z' G Ar^ 

Under the previous assumptions, there exists C > 0 such that the following bound 
holds for = Sq: 

(2.8) 'd^{n,K) > C\Bk\K^^ for all K G Z+, 

\n\>RKC 

where \ ■ \ denotes normalized Lebesgue measure on 9D. Thus, for eaeh p > 0, there 
exists a constant Cp > 0 such that 

(2.9) (K) > Cp\Bk\K^ for all K G Z+. 

Remark 2.2. Let us make a few remarks about the statement of Theorem O 
( 1 ) The assumptions of the theorem are translation-invariant. More precisely, 
suppose {an}n 6 Z is such that there exist C, 7 , and An as in the statement 
of the theorem. Then, if a'^ = for some k and all n, then there exists 
C' > 0 such that (12.71) holds (with the same choices for 7 and An)- This 
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is quite reasonable on physical grounds: global dynamical characteristics 
should not depend on how one chooses the origin. 

(2) The only step in the proof of Theorem l2. II which uses (12.61) is the step which 
passes between the Szegd and Gesztesy-Zinchenko cocycles (the GZ cocycle 
will be defined below). In particular, if a is not bounded away from 9D, 
then Theorem 12.II still holds, but with T replaced by Z. 

(3) Almost no modification is required to state and prove Theorem 12.11 in the 
half-line setting. In particular, one need only replace dSIZl) by 

sup max \\T{n,m-, z)\\ < CB?. 

z^Ar 0 <n,m<R 

The conclusions of Theorem 12.11 then hold verbatim in the half-line case. 
The modifications to the proof are straightforward and left to the reader. 

Notice that Theorem 12.11 supplies nontrivial dynamical bounds even when An 
consists of a single point for large R. Indeed, this gives a simple direct proof of 
time-averaged transport which is nearly ballistic (in the large p limit) for quan¬ 
tum walks in a Thue-Morse quasicrystal environment, and for quantum walks in a 
random polymer with critical spectral parameters, which we will describe in later 
sections. The latter example also shows that Theorem 12.11 sometimes can give 
much more information than results which relate dynamics to fractal dimensions 
of spectral measures, such as m Proposition 3.9]. Specifically, in random polymer 
quantum walks, one expects spectral measures to be of pure point type (it should 
be noted that no rigorous proof of this fact exists). Thus, the results of [11] give 
no information, while Theorem 12.11 yields almost-ballistic transport (at least in the 
limit p —> c»). 

2.3. The Gesztesy Zinchenko Cocycle and Dynamical Upper Bounds. In 

order to relate the Szegd cocycle to quantum dynamical transport behavior, we often 
need to relate the Szegd cocycle to solutions to the difference equation £u = zu 
with z S C \ {0} and u G C*. To that end, consider the matrices 

(2.10) P{a, z)=^- , Q(a, z) = ^ , a G D, z G C \ {0}, 

1 /2 

where p = (l — |ap) as before. Notice that 

det(P(a, z)) = det(Q(a, z)) = —1 for all a G D, z G C \ {0}. 

These matrices come from [3^ , though we follow the convention of [35] and replace 
a by —a in Gesztesy-Zinchenko’s definition (the definition of a GMV matrix in 
[25] also replaces our a by —a). One may use P and Q to capture the recursion 
described by the difference equation £u = zu m. the following sense. Any extended 
GMV matrix £ enjoys a factorization £ = CM. [43l Proposition 4.2.4], where C and 
M can be decomposed as direct sums of 2 x 2 unitary matrices of the form 



Specifically, 

£ = ^0(a2j), M = ^Q{a2j-i), 
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where ©(an) acts on coordinates n — 1 and n. If u G C* is such that Eu = zu, 
then we define v := C~^u. It follows that M.u = zv^ and we have E^v = zv. With 
$(n) := {u{n),v{n)Y, Gesztesy and Zinchenko show that 


( 2 . 11 ) 


^(n) 


P{an,z)^{n — 1 ) n is odd 
( 5 (a„, z)^{n — 1 ) n is even 


for all n G Z; see also m Proposition 3] for a generalization of this formalism 
to scattering zippers. This motivates the following definition. Denote Y{n, z) = 
P{an, z) when n is odd and Y{n, z) = Q(an, z) when n is even; then, the Gesztesy- 
Zinchenko cocycle is defined by 


{ Y{n — 1, z) ■ ■-Y^m, z) n>m 

I n = m 

Y{n, z)~^ ■ ■ ■Y{m — 1, z)~^ n<m 

If u, V, and $ are as above, we have 


(2.13) $(n — 1) = 2’(n, m; 2 ;)<I>(m — 1) for all n, m G Z 

by (12.111) : compare [25l Lemma 2.2]. We will often abbreviate the names and refer 
to Z{n, m; z) as a GZ matrix. 

By using the connection with solutions of the difference equation, (and hence 
with Green’s functions), we can use the GZ matrices to bound the tails of the wave 
packet; compare [l8l Theorem 7]. 

Theorem 2.3. Let E = Ea be an extended CMV matrix, consider the initial state 
i/) = ( 5 - 1 , and define 

Pr{N,K)= ^5^(n,iL), 

n>N 

Pi{N,K)= Y, a^{n,K). 

n< — N 


If a is bounded away from clD, there exists a constant Cq > 0 such that 

n2.Tr / o\ —1 

(2.14) Pr{N,K)<CoK^ 


max 

l<n<N 


(2.15) 


Pi{N,K)<CoK^ 


max 

l<-n<N 




dO 

2 \ de 

27r 


for all N, K G Z_|_. 


We use Ip = 6-1 in Theorem 12.31 purely for convenience; the result also holds for 
any finitely supported initial state (the constant Cq will depend on the choice of 
the initial state). 

We can also use the GZ matrices to bound the tails without averaging, at the 
cost of an additional power of k. 


Theorem 2.4. Let E = Ea be an extended CMV matrix, consider the initial state 
ip = S-i, and define 

PriN,K)= Y 

n>N 


(2.16) 
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(2.17) Pi{N,K)= Y, a^{n,K). 

n< — N 

If a is bounded away from 9D, there exists a eonstant Cq > 0 such that the left and 
right probabilities ean he hounded as follows: 

p'Ztz / . . o \ — 1 

(2.18) Pr{N,k)<Cok‘^ 

Jo 


max 

0 <n<N 


(2.19) 


Pi{N,k) < Cok^ 


max 

-N<n <0 


;,(„.0;e‘/-«)||y I 


2 \-^ dO 

2ti 


for all N,k £ Z+. 


These bounds on the tails imply corresponding bounds on the transport expo¬ 
nents in terms of bounds on the GZ matrices (averaged over the spectral parameter). 


Theorem 2.5. Given E = with a bounded away from 9D, suppose that there 
exist C £ (0, oo) and 7 £ (0,1) such that 




/.27r / 

(2.20) 

sup iF™ 

/ I max 

Kei,+ 

Jo \0<n<CKy 

and 


f2TV / 

(2.21) 

sup 

/ I max 

if ez+ 

/o \l<-n<Cif"' 


Z (n, 0; 


2 \ ^ de 

) 27r ^ 


2 \”^ de 

) ^ ^ 


for all m > 1. Then /3^(p) < 7 for every p > 0 and every finitely supported 

if £ eifL). 


Let us say a few words to compare and contrast Theorem 12.II with Theorem l2.5l 
First, there is a (very) rough physical heuristic which identifies the sizes of the 
transfer matrices with the sizes of the barriers that a wavepacket must tunnel 
through to escape compact regions. Thus, if all transfer matrices are sufficiently 
small, a wavepacket only observes small barriers, and hence, one expects it to tunnel 
effectively through all barriers; this vague statement is exactly what is made precise 
in the statement of Theorem nn In contrast, if transfer matrices are large, then 
a wavepacket observes large barriers and has much more difficulty tunneling out of 
compact regions; this statement is made precise in Theorem 12.51 

Moreover, notice that in contrast to Theorem l2.ll the application of Theorem [23] 
requires one to prove effective lower bounds on the transfer matrices on a set with 
nonzero Lebesgue measure. This contrasts with Theorem 12.II since as noted before, 
if one can prove effective upper bounds on transfer matrices at a single complex 
spectral parameter. Theorem 12. 1 1 already yields nontrivial conclusions! On physical 
grounds, this asymmetry in difficulty is unsurprising: to establish lower bounds on 
wavepacket propagation one need only prove that some portion of the wavepacket 
moves quickly, whereas to prove upper bounds one must simultaneously obtain 
effective control of the entire wavepacket. 


Remark 2.6. Let us conclude with remarks on a few extensions and variations on 
Theorems 12.31133 and 12.51 

(1) By general principles, /3j(p) < /3^(p) for any ip £ S and any p > 0. 
Specifically, if 7 > /3j(p), then |X|Ufc) ^ and it is not hard to use 
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this to show (K) < K'^P; compare the proof of [Ml Lemma 7.2], 

for example. In particular, the previous theorem immediately yields upper 
bounds on P^ip) for finitely supported '0 € £^(Z). 

(2) Over the course of the proofs, we use the hypothesis that a is bounded away 
from 915 in two places: to pass between Z(iV, 0;z) and 2'(2[iV/2j, 0; z) in 
the proof of Lemma [431 and, similarly, to pass between Z{N,Q]z) and 
2'(4[7V/4J, 0; z) in the proof of Lemma 14.91 Thus, one can modify the 
statements of Theorems mm and 12.51 to cover the case of general a’s 
by restricting n to be a multiple of 4 in the maxes which appear in Theo¬ 
rems 

(3) In the half-line case, the “left probabilities” make no sense, since a state 
cannot propogate to the left of the origin. However, the bounds on the 
right probabilities (equations (12.161) and (I2.18|) ') still hold true in the half¬ 
line setting. In particular, the correct generalization of Theorem [53] to the 
half-line setting is the following: if ( 12 . 201 ) holds (for some C, some 7 , and 
then all m > 1 ), then /3j(p) < 7 for all finitely supported ijj and all p > 0 . 

2.4. Motivation: Quantum Walks on Z. We now precisely describe quantum 
walks on the integer lattice and their relationship with extended CMV matrices, 
following [H]; see also [ 6 ]. A quantum walk is described by a unitary operator on 
the Hilbert space TL = .£^(Z)(g)C^, which models a state space in which a wave packet 
comes equipped with a “spin” at each integer site. As noted before, the elementary 
tensors of the form ( 8 ) e-|-, and Sn 0 64 , with n € Z comprise an orthonormal basis 
oi'H. A time-homogeneous quantum walk scenario is given as soon as coins 

(2.22) ®52)eU(2), nez, 

\Hn Hn / 

are specified. To avoid degenerate decoupling situations, we will always assume 
that ^ 0. As one passes from time t to time t + 1, the update rule of the 

quantum walk is as follows, 

(2.23) (5„ 0 e-i- qi^Sn+i 0 et -h ql^Sn-i 0 64 ,, 

(2.24) 0 64 , I—>■ ql^Sn+i 0 e-|- -|- q^Sn-i 0 64 ,. 

If we extend this by linearity and continuity to general elements of Ji, this defines a 
unitary operator U on H. Next, order the basis of H as in (EB , i.e. P 2 m — 

P 2 m+i = 0 64 , for m £ Z. In this ordered basis, the matrix representation of 

U :'H ^ H is given by 


•. 

0 0 (721 








ql \ 0 

0 







0 

0 


qf 






qi >^ 

0 

0 







0 

0 

21 

92 

q ¥ 






gP 

0 

0 


V 







•7 


which can be checked readily using the update rule (I2.23|) - (I2.24I) : compare [6l Sec¬ 
tion 4]. 
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We can connect quantum walks to CMV matrices using the following observation. 
If all Verblunsky coefficients with even index vanish, the extended CMV matrix in 
( 1 ^ becomes 


\ 


P 3 

0 

0 OiE Pb 
-03 0 0 

V ■■■ 

The matrix in (12.261) strongly resembles the matrix representation of U in (12.251) . 
Note, however, that > 0 for all n, so (12.251) and (j2.26l) may not match exactly 
when is not real and positive. However, this can be easily resolved by conju¬ 
gation with a suitable diagonal unitary, as shown in . Concretely, given U as in 
(I2.25L write 

fc e {1, 2}, r„ > 0, |a;^| = 1, 


(2.26) 


£ = 


0 

P-i 


0 

-Qf-l 


ai 

0 

0 

Pi 


Pi 

0 

0 

-ai 


03 

0 

0 

P 3 


define {A„}„gz by 

Ao = 1 , A_1 = 1 , X2n+2 = W^A2n, A2n+1 = W^A2n-l, 

and let A = diag(..., A_i, Ao, Ai,...). More precisely A : £^(Z) —>• £‘^{'1,) 0 maps 
Sn to Xn<Pn for each n € Z (where (pn is as in (I2.ip l. After a short calculation, we 
see that 

£ = A*C/A, 

where £ is the extended CMV matrix corresponding to the Verblunsky coefficients 


„ A 2 n -1 —^ A 2 n -|-1 12 ^ 

a2n = 0, Q; 2 n -|-1 = —;- Qn = “T- Qn > n G Z, 


( 2 - 27 ) — '^1 '^zn+l — , 'in — ^ In i 

^2n '^2n+2 

where equivalence of the two expressions for a2n+i follows from unitarity of Qn- 
Notice that the hypotheses ,q‘^^0 imply > 0 for all n G Z. 


The reader should be warned that many authors (e.g. [iHI]) use the transpose 
of the matrix U in (12.251) . However, since we imagine a vector i/; as a column 
vector with U acting on the left, our representation is more appropriate for the 
present paper. For quantum dynamical purposes, this is not a major concern for 
the following reason. If £ is an extended CMV matrix, then £^ = S*£'S, where S 
denotes the left shift on £^(Z) and £' has Verblunsky coefficients a'^ = an+i- Since 
all of our theorems are stable under shifting the Verblunsky coefficients, it does not 
matter whether one considers U or in what follows. 


Let us conclude with a few words on half-line quantum walks. A quantum walk 
on Z_|_ has state space %+ = (£^(Z+) 0 C^) 0 (do 0 64 .), and is specified by coins 
Qn £ U(2) for n > I. Additionally, one specifies a reflecting boundary condition at 
the origin, e.g., 


0 1 
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With this choice of boundary condition, one can check that (I2.23|) and (12.2411 define 
a unitary operator on and that this update rule has the form of a half-line 

CMV matrix with respect to the ordered basis {ipn ■ n € Z+j, where ipn is as in 

m- 


3. Power-Law Bounds on Szego Cocycles 


In this section, we will prove Theorem 12.II We begin by collecting a few useful 
facts. The following variant of the Parseval formula allows us to connect time- 
averages of observables to averages of the resolvent over the spectral parameter. 


Lemma 3.1. Suppose U is a unitary operator on Ji, if) G and {v?n} is an 
orthonormal basis ofT-L. For all K € and n, we have 


fc >0 


Jo 




de 

27r 


Proof. This is m Lemma 3.16]. 


□ 


Next, we will need a Gronwall-type variational estimate for the Szego matrices. 
The key ingredient in this estimate is furnished by the following useful variational 
lemma. Throughout the remainder of this section, let a be a fixed sequence of 
Verblunsky coefficients which is bounded away from i9ID). 

Lemma 3.2. Let z € (?D and w G C. Then, for n> 0, we have 

n—1 

(3.1) T{n,0-,w) = T{n,0-z)-Y^ {l - z-'^w) T{n,F, z)PT{£,0;w), 

£=0 

where 



Proof. This is [HI Proposition 1O.8.1J0 □ 

The following bound is an OPUC analog of [HI Lemma 2.1], which is the corre¬ 
sponding result for Schrodinger cocycles. 

Lemma 3.3. Suppose z G dlD and N > 1. Denote 

L{N) = max{jjT(n, m; z)]] : Jnj, jmj < N}. 

Then, for S G C with jdj < 1, we have 

(3.2) ||r (n, 0; ze^) || < L{N) exp {2L(N) ■ \n\ ■ J^J) 

for all n with |7 t,| < N. 


^In particular, the choice of Qq implies that (5o (g) is mapped to 5\ 0 e-j- so that U maps 'H+ 
to itself. 

^Actually, |44l Proposition 10.8.1] makes this statement only for u; € B, but the proof given 
there in fact establishes it for all € C. Alternatively, once one has (13.111 for all it; € B, the result 
for arbitrary it; € C follows immediately from analyticity. 
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Proof. Let us consider the case 0 < n < iV; the proof for —N < n < 0 is similar. 
By iU, we have 


n—1 

T (n, 0; ze^) = T{n, 0; z) - ^ (l - T(n, £; z)PT {(, 0; ze^) . 

£=0 

Now apply (ED to each T (£, 0; ze^) and iterate this procedure until all matrices 
correspond to spectral parameter z. Then, we obtain the estimate 


|T(n,0;.eO||<EQW+^ 




= L(7V)(l + 2L(iV)|(5|)'* 

< L{N) exp (2L{N)n\d\), 

where we have used ||P|| = 1 and that |1 — e'^j < 2\5\ for d S ID. 


□ 


With these tools, Theorem l2.1l can be deduced along the lines of [T7]. As noted in 
m, the Gesztesy-Zinchenko matrices are related to the Szegd matrices in a simple 
fashion. One has 


(3.3) z ^S{a,z)S{l3,z) = D^^P{a,z)Q{l3,z)Dz 

for all a, /? G ID and z G C \ {0}, where 


We warn the reader: this formula does not exactly match the corresponding identity 
in [12] ; this happens because [12] follows the conventions of [43] [44] with regards 
to the locations of the matrix elements of £, while we follow the conventions in 
[25] . Additionally, we take v = C~^u here, while m takes V = Aiu. With this 
relationship in hand, we can easily relate norm estimates on the two cocycles. 


Proof of Theorem \2.1\ Since M = sup„ |a„| < 1, we see that 
(3.4) 


Suppn < (1 — M'^) < OO. 


Using this and (13.311 . we see that there is a constant Co > 0 (which depends only 
on M = ||a||oo) such that 




for all n, m G Z, all K G and every z G 91D. Fix K G Z+, and let z G Bk be 
given. By definition, there exists z' G with \z—z'\ < l/K, and, by assumption, 
we have 

L = L{Rk, z') := sup ||T(n, m; z')|| < Ci?^ 

\n\,\m\<RK 

By Lemma EH we have the following for each n G Z such that 1 < n < Rk- 

2 


(u,0;ze^/*') 


< L exp \2L-n 


K 


< CR]^ exp 

















16 


DAVID DAMANIK, JAKE FILLMAN, AND DARREN C. ONG 


where we have used the power-law bound (j2.7|) in the second line. Since 
this yields 




< DR]^ for all 1 < n < Rk, 


where D = C'oC'e^'^. By (13.41) . we may enlarge I? in a {K, z)-independent fashion 
to get 

(3.5) Z < DR]^ for all 1 < n < Rk, 0 < j < 3. 

Now, take ijj = Sq, and define a pair of sequences u and vhy u := {£ — ze^^^) ip 
and V := C~^u; put $(n) = {u{n),v{n))^, as before. Throughout the remainder 
of the argument Ci, C 2 , ■ ■ ■ will denote constants which depend solely on £ (and 
neither on K nor z). Using (12.131) and (13.5L we have 

2 

(3.6) ||$(n)f > ^ ||$(j)f for all 1 < n < Rk- 

i=-i 

Notice that this uses det(2') = ±1 to get ||.^“^|| = \\Z\\ by the singular value 
decomposition. Since (£ — ze^/^)u = ip, we have 

a^p-iu-i - a^Pa-iUo + a^poui -I- PiPqU 2 - ze^^^uo = ipo = l- 
Consequently, we have 

2 

E > ^2. 

i=-i 

Summing (13.61) over ^RkI'^ < n < Rk, we obtain the lower bound 

E E mn)r > 

n>2RK /3 RK>n>2RK/3 

Since v = C~^u, we may use the explicit form of C to see that 

1 

|i'(?T-)P < ^4 E |m(?t-+ j)P for all n £ Z. 
i=-i 

Using this together with the definition of u, we obtain the following lower bound: 

2 


(3.7) 

n>RK /3 

By Lemma EH we have 


E (5n,{£ - = E 


R\-^^ 


n>RK /3 


2 

ap{n,K) = ^E 


^-2k/K^ 


'Ip {Rj k'j 


k>0 




de 

27r 


To complete the proof, denote = {9 G [0,27r) : e*® £ Bk}, and sum over 
n > Rk/3 to get 

2 


n>RK /3 


fs. ^ 


n>RK /3 



dO 


















SPREADING ESTIMATES FOR ID QUANTUM WALKS 


17 


which proves (12.81) . In the calculation above, we have used (13.71) to obtain the second 
line and K = to get the third. Applying (12.4|) to (12.8|) yields ()2.9|) . □ 


4. Upper Bounds on Spreading 


In this section, we will prove the bounds on left and right time-averaged prob¬ 
abilities from Theorem 12.31 The key ingredient is again Lemma 13.11 which gives 
us the following relation between each time-averaged single-site probability and a 
corresponding average of a matrix element of the resolvent off of cUD: 


(4.1) 


05.1 (n, AT) = —e 


2/K 


(Sn, 


E-e 


i/K+ie j 


5-1 


dte 

2t: 


Let’s take e = l/AT, and define G{z) = {£ — z) ^ for z G C \ a{£). Using (14.11) and 
the inequality < Stt, we then have 


(4.2) 

(4.3) 


Pr{N,K) < 3£ 
Pi{N,K) < 3 £ 



de, 

de, 


where we have defined 


(4.4) M,(iV,z)= ^ |(5„,G(z)5_i)p, 

n>N 

(4.5) Mi{N,z)= ^ |(5„,G(z)5_i)p. 

n< — N 


By symmetry, we will only prove the bounds on P ^, since the argument which proves 
the bounds on Pi is nearly identical. The overall strategy is to approximate £ by 
suitable sequences of “truncations” which converge strongly to £. One can then 
prove appropriate bounds on the quantities in (14.41) corresponding to the truncations 
iLemmas 14.21 and 14.3|) . Next, one controls the error incurred in passing between £ 
and the truncations ILemma 14.4p . Finally, one puts all of this together to bound 
Mr (and hence Pr) in terms of the GZ cocycle lLemma l4.5[l . The remainder of the 
section makes this discussion precise. 


Given iV G Z+, define four “truncated” CMV matrices £ 2 N, £ 2 N, E'an^ 
by choosing Verblunsky coefficients as follows: 


Ot2N{n) 


a2N{n) 



a(n) 

0 


3 

4 


when n < 2iV, 

when n > 2N and n even, 

when n > 2N and n odd, 

when n < 2N, 

when n > 2N and n even, 

when n > 2N and n odd, 
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{ a{n) 

0 

(_l)("+l)/23 

{ a{n) 

0 

(_l)(n-l)/23 


when n < 47V, 

when n > 47V and n even, 

when n > 47V and n odd, 

when n < 47V, 

when n > 47V and n even, 

when n > 47V and n odd. 


For convenience, let us also define 


.a = log(^). 

7?. = {z = £ C : 0 € [—7r/4,7r/4] U [37r/4, 57r/4], 0 < e < eo} , 

TZ' = {z = e C : 0 e [7r/4, 37r/4] U [57r/4, 77r/4], 0 < e: < Eq} ■ 


Let (p+(z, n), r+(z, n))^ be a solution to the GZ recursion subject to the initial 
conditions p+(z,—1) = z, r+(z,—1) = 1. Similarly, let (q+(z,n), S-i-(z,n))^ be a 
solution to the GZ recursion with q+(z, —1) = z, s+(z, —1) = —1. More precisely. 


(4.6) 


P+( 2 ,n) q+(z,n) 

r+(z,n) s+(z,n) 


= Z(n + 1, 0; z) 


z z 

1 -1 


where Z is as in (12.121) . Since det(P(a, z)) = det(Q(a, z)) = —1, we easily compute 
the Wronskian of these solutions: 


(4.7) W(z,n) :=p+(z,n)s+(z,n) - q+(z,n)r+(z,n) = 2(-l)”z. 


In particular, |IF| = 2|z| does not depend on n; moreover W is uniformly bounded 
on 77. U 77'. 

Let M±{z) denote the half-lattice Weyl-Titchmarsh m-functions corresponding 
to our operator £. Readers unfamiliar with the connection between these functions 
and the spectral theory of £ should consult Appendix A of [5S] . 

The functions M± are defined precisely in (2.137)-(2.140)], but here we 
will merely summarize the properties of M± most relevant to our proof. is a 
Garatheodory function, that is, an analytic function from ID) to the open right half 
plane. 71L_ is an anti-Caratheodory function, that is, an analytic function from D 
to the open left half plane. One of the reasons that these functions are important 
with regard to understanding the spectral theory of £ is that for z £ i9ID) we have 
this useful Herglotz representation, 

(4.8) M+{z)= f ^^^c7/r+(r), 

JdBT - 2 

where denotes the spectral measures of the restriction of £ to the right half-line 
(see [13 (2.123) and (2.137)]). There is also a similar formula for M_{z), but it 
is not relevant for our purposes. Using [13 (2.57), (2.61), (2.62), (3.5), and (3.7)], 
(splitting the GMV operator the way we do means setting /cq = — 1) we have 

(4.9) {{£ - z)~^S-i,6n) = d(z)s+(z, n) -I- 5(z)r+(z, n). 


for n > — 1, where 

(4.10) d{z) = 


1 + M_{z) 

2 z{M+{z) - M_{z)y 


b{z) 


{1 + M_{z))M+{z) 
2 z{M+{z) - M_{z)y 
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Let £ > 0 and z € C with log \z\ = £ be given. From (I4.6|l and (14.9p . we obtain 

(4.11) \b{z) - d{z)\ = (f - ^ < e-\ 

— 1 

The penultimate inequality is a standard consequence of the spectral theorem, while 
the final inequality follows from £ > 0. Notice that Re(M_|_(z)) < 0. Indeed, it is 
easy to see that 


Re 


T + Z 


T — Z 


l-\z[ 


T - Z 


< 0 


for every t with |r| = 1, so, by (ITSI) . we have Re(M+( 2 ;)) < 0. Consequently, 

'b{zy 


(4.12) 


arg 


d{z) 


= arg(M+(^))e (|,^ 


Combining (I4.11|) and (14.121) . we obtain 

(4.13) |6(z)| < £-^ and |d(z)| < £-\ 

Notice that the bound in (14.131) is a purely deterministic result. It holds for any 
choice of a and then for any z € C with log \z\ = e. 


Let us consider the GZ transfer matrices Y 2 N{'n,z) corresponding to the trun¬ 
cated Verblunsky coefficients a 2 N- Note that for n> N, we have 

>2^(211-b l,z)l2Ar(2n,z) = 

This matrix has determinant 1. We can calculate that the eigenvalues are 



(4.14) 


A± = A±(z) := 


2(z2 -b 1) ± V4 z 4 -b z2 + 4 

V7z 


We will wait until after Lemma lO to choose a branch of the square root 
yiz^ + z'^ + 4:. Notice that A± in our notation means something different than 
A± in [18]. After a short calculation, we see that the eigenvectors of 
l 2 Ar( 2 n -b 1, z)F 2 A( 2 n, z) are given by 


(4.15) 


= e^(z) = 


e?(^)\ f V7X±iz) - 4z 
ef{z)) 3 I 3 


for all n > N. 


Lemma 4.1. For all z € TZ, one of X±(z) has modulus strictly greater than 1, and 
the other has modulus strictly less than 1. In fact, A_ is bounded away from clD on 
TZ, i.e., there exists a constant cq < 1 with |A_(z)| < cq for all z € TZ. 

Proof. Consider £-oo, he., the 2-periodic extended CMV matrix with alternating 

Verblunsky coefficients ..., 0, —3/4,0, —3/4,_ Using standard calculations from 

Floquet theory, one can verify that 

E := ^(f-oo) = |z e c®> : |lm(z)| > . 

On the other hand, it is easy to see that z is a generalized eigenvalue of £-oo if and 
only if |A+(z)| = |A_(z)| = 1 (see [121 Section 3] for definitions and details). Thus, 
since 7^ (~l S = 0, the hrst statement of the lemma follows. To verify the second 
statement of the lemma, simply note that A-|_ and A_ are continuous functions of 
the spectral parameter away from z = 0 and that TZ is bounded away from E. □ 
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For iV G Z+ and \z\ > 1, define G 2 n{z) = {£ 2 N — z) 

Lemma 4.2. For n > 2N > 4 and z = G TZ, 

(4.16) \{ 6 ^,G 2 n{z)S.^)\<s-^ on liF 

je^ {z)p+{z, 2N -1)- T^[z, 2iV - 1)| 

(4.17) |(j.,g,^(z)j_i)|<£-^ on nr 

|ei {z)q+{z,2N - 1) - s+(z,2iV- 1)| 

Proof. Fix N G Z+ and consider a GZ recurrence corresponding to a 2 N- Let 
u = G2n{z)S-i = {£2n—z)~^S-i; thenu G £^(Z) and m satisfies = zu+S-i; in 
particular, u satisfies £ 2 nu = zu away from site —1. Let v = C~^u and 4) = (u, v)^, 
as usual (notice that v ^ P' also). We can write 

(4.18) 4-(27V-1) = (7+6++ C_e-, 

for some coefficients G± = C±{z,2N — 1) G C, where are defined in (j4.15|l . 
Using (12.131) . and 2N > 4, we have 

(4.19) 4>(2to - 1) = G+yff-^e+ + C.A^-^e" for all m>N. 

Note that (I4.19p doesn’t hold for $(2771), since are the eigenvalues of the two-step 
matrices Y 2 Ni‘ 2 n + 1, z)l2N(2n, z). However, using 

l2Ar(2TO, z) = for m>N, 

we must have u{2m) = v{2m — 1) and v{2m) = u{2m — 1) whenever m > N. 

Consequently, since the entries of 4> are in we must have C+ = 0. We thus 
have 

{ 62 n-i,G 2 n{z)S-i) = u{2N - 1) = C_(z,27V- l)ef{z), 

and 

{S 2 N: G 2 n{z) 6 -i) = u{2N) = v{2N - 1) = C_(z, 2N - l)ef{z). 

Using (SSI), we have 

(4.20) d 2 N{z)s+{z, 2N-1) + b 2 N{z)r+{z, 2N - 1) = C-{z, 2N - l)ef{z), 

(4.21) d 2 N{z)q+{z, 2N-l) + b 2 N{z)p+{z, 2N - 1) = C.(z, 2N - l)ef{z). 

where b 2 N and d 2 N are defined by (14.101) with £ replaced by £ 2 N- More specif¬ 
ically, notice that Z{n+ 1,0; z) only depends on a^,... ^an- In particular, the 
P+, ( 7 +, r_|_, s+-functions associated to £ and £ 2 n agree for — 1 < n < 2N — 1, which 
suffices to prove (14.201) by applying (14.91) . For (14.211) . notice that a 2 N{ 2 N) = 0 
implies 

sl^{z, 2N) = ql^ (z, 21V - 1) = g+(z, 2N - 1) 

rl^ (z, 21V) = p^'^(z, 21V - 1) = p+(z, 21V - 1), 

where the superscript 2N’s denote quantities associated to £ 2 n- Recall (14.7() : the 
Wronskian is given by 


= G TZ, 

\z\\X-{z)\^-^ 

er {z)p+izy 21V - 1) - r-+(z, 21V - 1) 
\z\\X-{zp-^ 

er iz)Q+iz, 2N -1) - s+(z, 21V - 1) 


W{z,n) :=p+{z,n)s+{z,n) - q+{z,n)r+{z,n) = 2(-l)”z. 
Solving (14.201) and (|4.21l) for C_, we get 


(4.22) 


W d2N 

efp+{2N-l)-efr+{2N-l) 
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(4.23) 


C_ = 


Wb 2 N 


—Cj q+(2N — 1) + 62 s+(2N — 1) 


In light of (I4.19p . this yields 


/A ^ \ 19 II Wd 2 NX^-^e^ 

{d 2 m-i,G 2 N[z)d_i} = u[ 2 m - 1) = — 


{S 2 m-l,G 2 N{z)S-i) = u{2m - 1) = 


erp+(2iV-l)-e2-r+(2IV-l) 
Wb 2 NX-~^e^ 


-e^q+(2N - 1) + e^s+(2fV - 1) 

The statement of the lemma for n > 2N odd then follows immediately from (I4T1) . 
(I4.13[) . and (14.1511 . The statement for even n > 2N follows from 

{b 2 m,G 2 N{z) 6 -i) = u{2m) = v{2m — 1) for all m> N. 

□ 

Using a completely analogous calculation, we obtain the following lemma for the 
Green’s function of £. 

.3. Define G 2 n{z) = {£ 2 N — z^ for \z\ > 1. Then, for z = G TZ 


Lemma 4 

and n > 2N > 4, 


(4.24) 

(4.25) 


G2N{z)b-i 
bn, G2 n{,z)5-1 


<£-'■ 


|z||A(^)|t 


-N 


Cl {z)p+{2N - I, z) + r+{2N - l,z)\ 
\z\\X-{z)\^-^ 


e^ {z)q+{2N -l,z) + s+(2N - 1,2:)| 

Given iV G let us denote N_ = 2[iV/2j. For z ^ 9D, dehne 

2 


S{N,z)= Y. |(,5„,Gjv(^)<5-i) 


n—N-\-l 
00 

S{N,z)= Y |(^n,Gw(z)<5_i 

n=N+l 

Lemma 4.4. For all z = G TZ, we have 

(4.26) e^S{N, z) < Mr{N, z) < e-'^S{N, z), 

(4.27) e'^S{N, z) < Mr{N, z) < e-'^S{N, z). 

Proof. Using a standard resolvent identity (e.g. [48l Lemma 6.5]), we get 
{£ - z)-i - (% - ^)-^ = -(f - z)-^ {£ - £n) {£n - z)-\ 
Since ||(f — 2:)“^|| < s~^ we then have 

00 

Mr{N,z)= Y l(^n,G(z)J_i)p 

n—N+1 

00 

= ^ \{5n,GN{z)5.i-G{z){£-£^GN{z)5-i)\’^ 

n=N+l 

00 

<2S{N,z) + 2e-m£-£N)f Y |('5n,Gw(z)<5_i)f 

n—N-\-l 
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<e-^S{N,z). 

In the first inequality, we have used |a + < 2|ap + for complex numbers 

a and 6 , which follows readily from Cauchy-Schwarz. The proof is similar for the 
other inequality. The proof in the S case is identical. □ 

Lemma 4.5. For all z = G TZ, we have 
(4.28) Mr{N,z)<e~^( max ||Z(n, 0; z)||^ 

\4:<n<N 

The implicit constants depend on M = ||q;||oo 7 but nothing else. 



Proof. The previous lemmas show us that 
(4.29) Mr{N,z)<e-‘^ 


Mp+(iV-l,z)-r+(iV-l, 2 :) 


(4.30) Mr{N,z)<e-'^ 


3 

^A_(z) - P+iN-l,z) + r+{N-l,z) 


-2 


-2 


-2 


(4.31) Mr{N,z)<e ^ ^^A_(z)-^z ^ j g+(iV - 1, z) - s+(iV - 1, z) 

(4.32) MriN,z) <e~^ ^^A_(z) - g+(iV - 1, z) + s+(iV - 1, 2 ) 

Notice that (I4.29|) and (I4.31|) implicitly use the claim from Lemma|4T]that | A_ | < cq 
on TZ to uniformly bound 


00 


(4.33) ^|A_(z)P<- -< 

7^0 1-^0 

for all z G TZ. This step explains the need for several truncations, since we cannot 
have a bound like (j4.33|) on TZ', which means that we lose control of the implicit 
constants in that region with this pair of truncations. Let us denote 

(y/T 4 ^ 

(4.34) U={ ' 


(4.35) 


U = 


^^A_(z)--z J p+(iV - l,z) - r+(7V-l,z), 
P+iN-l,z) - r+{N-l,z). 


Adding and subtracting the two equations we note that 

\u\ + \U\>1 \VlX-{z) - 4z-i| b+(7V - l,z)|, 

and also 

|[/| + |t/|>2|r+(iV-l,z)|. 

Since |A_| < 1 and 1 < |z| < we must have •\/7A_(z) — 4z“^ 7 ^ 0. 

This implies 


(4.36) M.(iV,z) < < £-4 (|p+(^- 1,^)|2 + 
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We can similarly also show that 

(4.37) Mr{N,z) < e-^ {\q+{N- l,z)|2 + |s+(iV - 1,. 

For the sake of notational convenience, in the following calculation let us use 
p+,q+,r+, s+ as a shorthand forp+(^—1, z), g+(fV—1, z), r+(^—1, z), s+(iV—1, z). 
We can use (14.61) to observe that 


\\Z{N,0;z)\\'^ = sup 

(A,B)/(0,0) 


< sup 
(A.B)/(0,0) 


A 


P+ 

r+ 


B 


B 


z 

-1 


A 

{\A\^ + \B\^){\p+\^ + \q+\^ + \r+\^ + \s+\^) 


\A + B\^ + \A-B\' 


= -(|p+P + |g+p + |,+ |2 + |,^|2) 

< max(|p+p + |r+p, \q+\^ + |s+p). 

Since a is bounded away from 5ID), we have Z{N, 0; z) < ^(fV, 0; z), which implies 

MriN,z)<e-^\\Z{N,0;z)\\-^. 

Observe that Mr{n, z) > Mr{N, z) whenever n < N. Thus, the previous argument 
gives the following for 4 < n < iV: 

MriN, z) < Mr{n, z) < £-4 {\\Z{n, 0; z)f)-' . 

Since the preceding inequality holds uniformly in 4 < n < iV, we get (I4.28|) . □ 


We want to extend all these results to arbitrary 0 G [0, 27r) by using the second 
pair of truncations of E. We can then obtain a new version of Lemmas 14.21 and 14.31 
Notice that the four-step GZ matrices of the truncated CMV operators take the 
form 


Z4^(4n-|-4,4n;z) 
^4Ar(4n-|-4,4n;z) 


11 

f 16z^ — 9 

12(z-z-i) 

7 ' 

l^-12(z- z-i) 

16z-2 - 9 

1 1 

f 16z^ — 9 

-12(z-z-i) 

7 * 

1^12(z- z-i) 

16z“^ — 9 


for n> N. Both have eigenvalues 


, , 8(z 2 -h z-2) - 9 ± ^/(8(z2-h z-2) - 9)2 - 49 

- 

The first has eigenvectors v± = given by 

- z“i) ~ ^ ^ 

The second has eigenvector (—v^, 1)^. Following previous arguments, we can prove 
bounds for the decay of the corresponding Green’s functions as before. 


Lemma 4.6. Define G'^j^{z) = — z) ^ for \z\ > 1. For n > AN > 8 and 

z = G 


(4.38) 


|((5„,g;^(z)<5_i)| <£ 


-1 


l^l|A'(^)|? 


-AT 


'(z)p+(z,4iV — 1) — r+(z,4iV — 1) 
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(4.39) 

|(<5„,GU(^)<5-i)| 

< 

kl|A'(^)l^ 

-N 

r 

h{2)Q+(.z,4:N- 1) - 

s+{z,4N- 1) ■ 

Lemma 4.7. Define G'^j^{z) 
z = G TZ!, 

= 

— z^ for \z\ > 1. 

For n > 4:N > 

(4.40) 

(5„,G'4^(z)5_i) 

< e-U 

|z||A'(^)l^-^ 


vf {z)p+{4:N - l,z) + r+{4:N - l,z) 

(4.41) 

(<5„,G'4^(z),5_i) 

< F-U 

kl|A'(^)|- 

i-N 


(z)g+(47V- l,z) + s+(47V- l,z) 

For z ^ 9D, define 





S'{N,z) = 

oo 

E 

(dn, G'4[N/4J (^)A-l^ 

2 

5 


S'{N,z) = 

E 

^i5„,G4|^^/4j(z)(5_i^ 

2 




A new version of Lemma 14.41 is also immediate. 

Lemma 4.8. For z = € TZ', we have 

(4.42) e^S'{N, z) < Mr{N, z) < e-'^S'{N, z), 

(4.43) e^S'{N, z) < Mr{N, z) < e-'^S'{N, z). 


Notice that Vi{z) is bounded away from zero for z G TZ'. To see this, notice that 
(0,1)^ is an eigenvector of Z'^j^{4:N + 4,47V; z) if and only if z — z~^ = 0, that is 
z = ±1. Since TZ' is bounded away from {±1}, it follows that u)“(z) is bounded 
away from zero. This observation allows us to extend the arguments (and hence 
the conclusions) of Lemma |43] to spectral parameters in TZ'. 

Lemma 4.9. For z = € TZ', we have 

(4.44) Mr{N, z) < ( max ||Z(n, 0; z)||^ 

\4<n<N 

The implicit constants depend on M = ||a||oo) but nothing else. 

Proof of Theorem \‘d.3i Follows immediately from Lemmas 14.51 and 14.91 except that 
we have the max in (12.141) running over n > 1 instead of n > 4. Since a is 
bounded away from 9D, making the appropriate adjustments only modifies the 
implicit constants in an iV-independent fashion, so this is hne. 

□ 



5. Upper Bounds Without Averaging 

By using techniques from complex analysis and a version of the residue calculus 
for normal operators, we can prove variants of the upper bounds from the previous 
section without averaging in time; compare |19j . 

Lemma 5.1. Let £ he an extended CMV matrix. We have 

(5.1) (5„,f^A_i) = -^^z'=(<5„,(f-z)-i5_i) dz 
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for every fc € Z, all n £ Z, and every positively oriented contour F m C which 
encloses D. 


Proof. This follows from the Dunford functional calculus for normal operators. See 

EaEiiini. □ 


Lemma 5.2. Let £ he an extended CMV matrix, and consider the unaveraged left 
and right probabilities Pr and Pi. We have 


(5.2) Pr{N,k)< f V 'd_i\ 

■^0 n>N \ ^ ^ ' 

(5.3) Pi{N,k)<j^J Y. {5r,,(£-e^/^+^^y^ 


d6 

27r 


dO 


for all k £ 'L \ {0}. 

Proof. Consider the contour 

r = {exp +it) : 0 <t < 27r} . 
We apply Lemma [5T] and Jensen’s inequality to get 
P^{N,k) = Y 


n>N 


= Y \{Sn,£^sY\' 

n>N 

E r 


n>N ‘ 




d0 


which proves (15.21) . The proof of (15.31) is nearly identical. 


□ 


With these pieces in place, we are now able to prove Theorem 12.41 very easily. 


Proof of Theorem \2.4\ Follow the proof of Theorem l2.31 but replace (14.21) and (14.31) 
with Lemma [5.21 More specfically, we have 


Pr{N,k) < 



Mr (iV, e^/^+y 


dO 

27r 


< M 


max 

0<n<Af 




dO 

2?! 


The first inequality is (15.21) . and the second is Lemma IT3] The bound on Pi works 
similarly. □ 


We can use this to prove the bound on the /3’s from Theorem 12.51 


Proof of Theorem \2.5\ For = (5_i, this follows immediately from Theorem 12.41 
More precisely, the assumption (I2.20|) together with Theorem 12.41 implies that 
Pr{Ck~^,k) + Pi{Ck~^,k) goes to zero as fc —>■ oo faster than for any s > 0. 
Enlarging C if necessary, we may assume that a.,p(n, k) = 0 whenever \n\ > Ck + C. 
This implies 

= E ( 1 ^ 1 ^ + 
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= {\n\^ + ^)ai;in,k) + ^ {\n\P + 1) a^{n,k) 

|n|<CfcT Ck'i <\n\<Ck+C 

< Ck'^P + l + {Ck + C)P {Pr{Ck"<, k) + Pi{Ck'^,k)). 

for all k G Z+ and p > 0. Consequently, we obtain the theorem for ip = (5_i. For 
Ip = Sn with n gIj, the theorem follows by applying the case ip = S_i to the CMV 
operator where S : £^(Z) —>• £^(Z) denotes the left shift. The result 

for finitely supported ip's is then obvious. 

□ 


Part 2. Examples 

6. Delocalization for Quantum Walks in Polymers 

6.1. Setting. We describe a quantum walk analog of polymer models; compare 
[28] . Fix a finite subset A C U(2), and let A* = denote the free monoid 

over .A, that is, the set of all finite words which can be obtained by concatenating 
elements of A. A polymer model is given as soon as one chooses n> 2 and elements 
ui,... ,Un G A*, which we call the basic chains of the model. We can then define 
the family of quantum walks generated by these basic chains to be those whose 
sequences of coins may be obtained by concatenating elements of {mi,...,u„}- 
More formally, let D = {1,..., n}^, and, for each w € D, let Uui be a quantum walk 
obtained by concatenating ... ,Ucj_^,Ucjg,Uuji,... with the position of the origin 
normalized to be at the start of For each j let Tj{z) denote the Szego transfer 
matrix across the coin sequence Uj at spectral parameter z. We say that a spectral 
parameter z G c® is critical for the random polymer model if the following three 
conditions hold: 

(1) |tr(Tj( 2 :))| < 2 for every 1 < j < n. 

(2) Tjlz)Tk{z) = Tk{z)Tj{z) for all 1 < j, fc < n. 

(3) If |tr(r/(z))| = 2, then Tj{z) = ±1. 

We distinguish such spectral parameters as “critical” in analogy with the critical 
energies for the Schrodinger equation with random polymer potential; in the context 
of the Schrodinger equation, these correspond to energies at which localization 
lengths diverge and Anderson Localization breaks down [17], [21], [^, [28]. At 
each of these critical values, it is easy to produce a uniform bound on the associated 
transfer matrices; by general arguments m, this ensures that the critical energy 
is in the spectrum. 

Similarly, in our CMV operator setting any critical spectral parameter is con¬ 
tained in the spectrum of every operator of the form with w S D. 

Example 6.1. It is not too hard to construct polymer models with critical spectral 
parameters. For example, notice that 

S{a, 1)5(0, l)5(-a, 1)5(0,1) = / 

for all a S ID). In particular, z = 1 is a critical spectral parameter for any random 
polymer model with basic chains of the form {Rg,R-g) with 9 G (—§, f), where 
Rg denotes the rotation matrix defined by 

Pg= 


cos 9 — sin 9 
sin 9 cos 9 
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In Section El we will study the a quantum walk modelled on the Thue-Morse 
subshift. This can also be incorporated into the framework of polymer models, and 
we will see that this model has many critical spectral parameters. 

Theorem 6.2. Given a polymer model with a critical spectral parameter, one has 
> 1 - 

Proof. Let zq denote a critical parameter, and put An = {zq} for all R > 0. Using 
the definition of criticality, we can choose Co > 0 so that (12.71) holds for all i? > 0 
with 7 = 0. Since \Bk\ = 2K~^, we deduce (K) > KP~^ from (12.9p . Thus, 

we get > 1 - i. □ 


7. The Fibonacci Quantum Walk 

7.1. Setting. In this section, we will discuss a time-homogeneous quantum walk 
on Z with two coins, whose distribution on the line is modulated by the Fibonacci 
subshift. We will apply both of our general theorems to deduce upper and lower 
bounds on the corresponding quantum walk. This is by far the most substantial 
application of our methods; it will require the most effort in order to verify the 
inputs needed to apply the general theorems. The key observation here is that 
one can restrict the range of parameters in such a way that the spectra of the 
canonical periodic approximants exhibit the same combinatorial structure that one 
sees in the Schrddinger case; compare [39] . Since our goal is to demonstrate the use 
of the transfer matrix method in quantum dynamics, we will freely make various 
simplifications throughout this section in order to increase clarity (without removing 
any of the essential obstacles). In particular, we do not attempt to locate optimal 
parameter ranges, nor do we attempt to find optimal constants in our proofs. 

First, we recall how the Fibonacci sequence and subshift are generated. Consider 
an alphabet A with two symbols, A = {a, 6}, and let A* denote the free monoid 
generated by A, i.e., the set of all finite words obtained by concatenating elements 
of A. The Fibonacci substitution S sends a to ab and b to a. Naturally, S can 
be extended by concatenation to maps A* —>■ A* and —>■ AA^ which we will 

also denote by S. The Fibonacci sequence, which we denote by up, is the unique 
element of A?"+ that is invariant under S. It is the limit of the words Sk '■= S^{a) 
as fc —oo in an obvious sense. More precisely, with sq = a, si = ab, S 2 = aba, 
etc., we see that Sk is a prefix of for all k, so we define up to be the unique 
sequence in A^+ which has Sfc as a prefix for every n > 0: 

Up = abaababaabaab .... 

The Fibonacci subshift Up is then defined to be the dynamical hull of up in A^, 
the space of two-sided sequences over M; more precisely. 

Up = {w G : every finite subword of ui occurs in up} . 

To generate a quantum walk in a Fibonacci environment, choose 9a,9b G (~§i §)j 
and consider the rotations 

f cos 9a -sin6(a\ f cos 9b -sm9b\ 

- [sm9a cos9a)^ [sin 9b cos 9b ) ' 

Given uj G Up, we associate a sequence of coins {Quj,n}ne 7 . via Quj,n = Ru^- The 
associated unitary operator will be denoted by Uui■ Inspecting (12.2711 one sees that 
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Uu> already has the form of an extended CMV matrix with Verblunsky coefficients 

Q: 2 n +1 = sin(0^^), a 2 n = 0, for all n £ Z. 

We will therefore denote Uuj by to emphasize this fact. Notice also that we 
must avoid ±7r/2, for if either 9a or 0i, is ±7r/2, then we have |a 2 n+i| = 1 for 
infinitely many n £ Z, so decouples into a direct sum of finite blocks, which 
means that the dynamics are trivially localized. By a standard argument using 
minimality and strong operator approximation, there is a uniform compact set 
E = J^{9a,9b) C 9D such that E = a{£ai) for every lo S flp- Specifically, if 
a;i,a ;2 £ ^f, then wi may be approximated by shifts of 012 in the product topology, 
and thus, may be approximated in the strong operator topology by operators 
which are unitarily equivalent to £ 012 - One may then apply standard results (e.g. 
[40l Theorem VIII.24]) to deduce a{£uji) ^ a{£ai 2 )', equality of spectra follows by 
symmetry (one may exchange the roles of wi and W 2 in the foregoing argument). 

To avoid cluttering the notation, we will suppress the dependence of the various 
objects in this section on the choice of 0a and 9b. 

We consider the initial state ip = 5^) and study the spreading in space of £Jp%p as 
\n\ —^ 00 with respect to this basis. For the sake of clarity and simplicity, we will 
focus on the dynamics generated by the element uiq £ obtained by applying 
iteratively to the germ h\a, where the vertical bar separates sites —1 and 0. We 
may write wq = Upa6|uF, where denotes the reversal of the word u. Thus, wq 
agrees with uf on the nonnegative semiaxis, and it has the reflection symmetry 

(7.1) uJo{—n) = uJo{n — 3) for all n € Z \ {1, 2}. 


Equivalently, one may write 


(7.2) 


a;o(n) 


a {{n + l)(p 1} e [l-i^ \1) 

b {(n + l)(^-i} e [0, ’ 


where {x} = x — [xj denotes the fractional part of x G R and denotes 

the golden ratio. One can deduce (m from (4] Lemma 1], for example. 

The structure of the remainder of the section is as follows: In Subsection 17.21 we 
describe the trace-map formalism for the operators £a ;, and we state our dynamical 
bounds precisely in terms of quantities associated to the trace map. Subsection 17.31 
contains a complete proof of the main theorem, modulo a pair of lemmas which con¬ 
tain some technical estimates on the trace map and on the growth of the Gesztesy- 
Zinchenko cocycle for this model. Subsection 17.41 works out a version of Raymond’s 
combinatorial analysis of the spectrum in this setting. Finally, in Subsection 17.51 
we prove the technical estimates from Lemmas 17.81 and Lemma 17.91 


7.2. Trace-Map Formalism and Pricke-Vogt Invariant. In this section, we 
fix a = a^jo and £q = £a^^ ■ The renormalization map for the Fibonacci sequence 
will play a key role in our analysis. More concretely, define 

M„(z) = Z(2F„, 0; z), n > 0, z G C \ {0}, 

where Z denotes the Gesztesy-Zinchenko cocycle as before, and Fn denotes the nth 
Fibonacci number, normalized by F-i = Fq = 1- We also define 
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The hierarchical structure of the Fibonacci substitution word produces a recursive 
relationship amongst these matrices, viz. 

(7.3) Mn+i{z) = Mn-i{z)Mn{z) for all n > 0 and all z € C. 

This recursion leads to a number of nice consequences, which are well-known and not 
hard to prove; for the corresponding statements and proofs in the context of the dis¬ 
crete Schrodinger equation, consult [47l Proposition l(ii)]. By the Cayley-Hamilton 
Theorem, (1731) leads to a recursion amongst the half-traces Xn ■= ^tr{Mn{z)), 
namely, 


(7.4) Xn +2 = 2XnXn+l — Xn -1 for all n > 0. 

As a consequence of this recursion, the traces have a first integral given by the 
so-called Fricke-Vogt invariant. More precisely, if we define 

I{u, V, w) = + v'^ + w'^ — 2uvw — 1, 


then 


(7.5) I{Xn+l,Xn,Xn-l) = I{Xn+2, Xn+1, Xn) for all n > 0. 

Consequently, by abusing notation, we will write I{z) := I{xk+i{z), Xk{z), Xk-i{z)) 
for k > 0 and z £ C. Evidently, Mk{z) is the monodromy matrix for a periodic 
operator whose sequence of coins is defined by the sequence ujk £ which repeats 
Sfc periodically. By Floquet theory, the spectra of these periodic operators are given 
by 

(7.6) (Tfe := {z e i9ID) : |a;fc(2:)| < 1}, k> -1. 

For more details on Floquet theory and the spectral characteristics of CMV matrices 
with periodic coefficients, the reader is referred to [44l Section 11.2]. In this setting, 
the role of the coupling constant may be played by 

(7.7) = fj.{9a,0b) ■■= inf min/(z). 

fc> —1 z£(Tk 

Let us also dehne 


(7.8) K = k(6(q, 0b) = I sec6a tan 0b — tan0a sec 0b| = | sec0a|| sec0b|| sin0a — sin0b|. 

We can now state our main theorem, which describes the dynamics defined by Sq 
in the quantum walk analog of the large coupling limit. 


Theorem 7.1. Let tt/A < 0b < 7r/2 be given. There exist constants m = m(0b), 
M = M{9b), and X = A(0b) such that if db < 0a < 7r/2 with fi > X, one has 


(7.9) 


1 


3t + r] p — 3 t — r] 


1-I-T p(I-I-t) p(l-l-r) 


< 0So(P) ^ I^SoiP) ^ 


21ogy> 
log^ ’ 


where 

(7.10) £, = my/JI, 


= r] 


login. ^ 2 log ((k-|-2)(2k-|-5)^) 

log ip ' log 


and K is as in (17.81) . In particular, {p) goes to zero at least as fast as 
constant/ log as p ^ oo. 


Remark 7.2. A few remarks are in order. 

(1) The conclusion of the theorem likely still holds with p = min^ eT,I{z). 
However, our choice of p makes the proofs somewhat easier, so we use dZlD 
to make the exposition more digestible. 
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(2) Our method of proof for the lower bounds can allow 0 < < 7r/4, but the 

proof of upper bounds breaks at Ob = tt/A. However, on heuristic grounds, 
one may still expect that /3 \ 0 as da ^ 7i/2 for other values of Ob- 

(3) For any fc > 0, recall that Sk = S^{a)-, additionally, define := S^{b). 
By [101 Lemma 3.2], every ui G Op may be decomposed as a concatenation 
of subfactors of the form Sk and s], in a unique way; this is known as the 
k-partition of ui. Theorem [Td] can be extended to arbitrary elements of flp 
by analyzing various possibilities for the fc-partition near the origin, but 
we will not present the details, since this clutters the presentation without 
adding substantial new content. 


7.3. Proof of Theorem 17.11 First, for Theorem 1 7. II to have any interesting con¬ 
tent, we need to verify that p, may indeed be made arbitrarily large by suitably 
choosing 9a and Ob- To that end, let us begin by computing / as a function of z. 
This is done in |16j . but there are a couple of typos and a sign error, and, most 
importantly, we have a substantially simpler expression for I{z\ For these reasons, 
we produce a corrected calculation here for the official recordO We have 


Mo{z) = sec 9a 


Ml {z) = sec 9a sec 9b 
Thus, for z £ 9D, we have 


z — sm b 

— sin da 

z^ -I- sin 9a sin Ob 


* 

sin da sin df, -|- z 



(7.11) x_i(z) = Re(z) secdf,, 3 : 0 ( 2 :) = Re(z) sec da, 

(7.12) 3 : 1 ( 2 ;) = Re(z^) sec da sec df, -I- tan da land;,. 

Consequently, 

/(z) = Re(z)^(sec^ da -|- sec^ 9b) + (Re(z^) sec da sec Ob + tan da tand;,)^ 
— 2Re(z)^ sec 0a sec 9b (Re(z^) sec da sec 9b + tan da tan 9b) — 1 
for all z G 9D. 


Lemma 7.3. For every z £ 9D, we have 
(7.13) /(z) = K^(lm(z))^. 

Proof. Let z = £ 9D be given. To compactify notation a bit, denote 

Ca = cosda, Pa = Secda, Sa = sillda, ta = tan da 

Cb = cos9b, pb = sec Ob, Sb = srnOb, = tandf, 
c = Re(z) = cost, s = lm(z) = sint, C 2 = Re(z^) = cos(2t). 

We have seen above that 

7(z) = C^{pl + pI) + {c 2 PaPb + tatb)^ “ pIpI (C 2 + SaSb) - 1. 

After liberally using trigonometric identities, one can deduce ()7.13l) . The detailed 
calculation follows: 

I = C^{pI + pI) + {c2PaPb + tatbf - pIpI (c 2 + SaSb) - 1 


^Note that m uses the Szego transfer matrices Mn{z) = z ^■^T(2Fn, 0; 2 ), and we use the 
GZ transfer matrices. In light of lira . this does not affect the value of Xk(z) for z € and 
therefore does not affect I{z) for such 2 :. 
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— PaPb ”1" + (C2 + SaSb)^ — 2c^C2 — 2c^SaSt,) — 1 

= pIpI (c^(Ca + cl) - C 2 + slsl - 2s^SaSb) - 1 
= pIpI ((1 - S^)(Ca + cl) + 2s^ - 1 + (1 - cl){l - cl) - 2s^SaSb) - 1 

= pIpI (s^(Sa + - 2s^SaSf, + C^C^) - 1 

= plpls^[s,-Sb)^ 

= (Im( 2 ))^ (sec Pa tan 0b — tan 0a sec0b)^ , 

which proves (17.131) . The third line is obtained by expanding the binomial square 
and using the identities C 2 = 2c^ — 1 = . The fourth line uses C 2 = 1 — 2s^ and 

the Pythagorean identity. The fifth line is simple algebra (and uses the Pythagorean 
identities some more). □ 

Remark 7.4. From the expression in (17.131) . we make several observations. 

(1) The invariant is always nonnegative on the spectrum. Indeed, for any 
z G 9ID), I(z) > 0. 

(2) One cannot make the invariant uniformly large on clD by suitably choosing 
0a, 0b- Indeed, one always has /(I) = /(—I) = 0. 

(3) For any fixed 0b, one has 

lim K(0a,0b)=OO. 
e„y7i-/2 

Thus, for fixed 0b, we would like to think of 0a ~ 7r/2 as the quantum walk 
analog of the large coupling regime. We will see later that one can force 
the spectrum away from ±1 for some parameters, which allows us to make 
the invariant uniformly large on S. 

(4) Similarly, we think of 0a ~ 0b as the “small coupling regime,” since we can 
make I uniformly small on 9D by making 0a and 0b sufficiently close. 

Let us note that we can characterize the uniform spectrum of Soj as the dynamical 
spectrum, i.e., the set of complex numbers at which the associated trace map has 
a bounded orbit. More precisely: 

Proposition 7.5. Fix 0a, 0b G (~f ^ f)^ letT, denote the uniform spectrum of every 
Eo) with uj G Hf, and put 

Boo = {z G C : (xn(z))^^_i is a bounded sequence} . 

Then S = Boo- 

Proof. The inclusion E C Boo is in m Lemma 2.3]. The other inclusion can be 
proved using ideas from lochum-Testard m- Namely, given z G Boo, one can 
linearize the matrix recursion as in m and use the assumption of bounded traces 
to derive a power-law upper bound on the norms of the matrices ||.^(n, 0; z)|| with 
n G Z. Indeed, this proof is nearly identical to the proof of Lemma 17.91 Once one 
has a power-law bound on the growth of the transfer matrices, the 2 in question 
belongs to the spectrum by general principles (e.g., [m Theorem 6]). □ 

In order to decide whether or not a given z G 9D is a member of Boo, it is helpful 
to have a computable condition which guarantees escape of the corresponding trace 
orbit. The following statement is not hard to prove; it is m Lemma 4.3] applied 
to the special case a = ip~^. 
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Lemma 7.6. Fix z G C and S > 0. If there exists fco > 0 such that 

(7.14) |a;fco(z)| > 1 + 5, |a;feo+i(z)| > 1 + 5, and \xk„+i{z)\ > \xkQ-i{,z)\, 

then the sequence {xk{z))'^_i is unbounded. Moreover, if (17.141) holds, then 
\xko+jiz)\ > (1 + 5/2)^^-! for all j > 0. 

We are now able to verify that the Fricke-Vogt invariant can be made uniformly 
large on the spectrum in the large coupling regime. 


Proposition 7.7. For any A > 0, there exist choices of 9a, 0b > 0 such that fJ.> X. 
More specifically, for any 0 < 9b < there exists (fo = 4>Q{9b) > 9b such that 
I{z) > A for all z € ak and all k > —1 whenever < 9a < ^. 

Proof. We will show that we can force the spectrum into the cone where the in¬ 
variant is large by forcing the escape condition to hold for ko = 0 and 5 = 0 in a 
suitable cone containing the real axis. Fix 0b > 0, and suppose 9a > 9b is large 
enough that 

( 7 . 15 ) = {sec 6a t&n 9b — tan 9a sec 6b)'^ > — ^ — 

sin 9b 

( 7 . 16 ) (2cos^0b — 1 ) + sin 0(1 sin > cos 0a- 


Notice that the second condition can be satisfied by taking 0a sufficiently close to 
+7r/2, since 

2 cos^ 06 — 1 + sin 06 > cos^ 06 > 0. 


Now, recall: 


x-i{z) = Re(z) sec06, xo{z) = Re( 2 ) sec0Q, 
xi{z) = Re(z^)sec0a sec06 + tan0a tan06. 

It is easy to see that 

\x-i{z)\, \xo{z)\ > 1 

whenever |Re( 2 :)| > cos06, since 9a > 9b > 0. We also want to show |a;i| > |a:_i|, 
which holds whenever 

(7.17) |Re(z^) + sin0Q sin06| > |Re( 2 ;)| cos 0 q. 

By (|7.16l) . we see that (I7.17|) holds for all 2 with |Re( 2 )| > cos06. Since |x_i( 2 ;)| > 1 
for such 2 , we also get |a;i(z)| > 1 for those values of z, so the escape condition 
(j7.14p holds for fco = 0, 5 = 0, and all z G i9ID) with |Re(z)| > cos06. It follows that 

(Tfe C a-I = {z € 9D : |Im(z)| > sin 06} for all k > —1. 

Consequently, we have I{z) > X on ak for all /c > —1 by (j7.13p and (17.151) . □ 

The proof of Theorem 17.11 relies on the following pair of technical lemmas. The 
first lemma establishes upper and lower bounds on the derivatives of the iterates 
of the trace map, while the second proves power-law bounds on the growth of the 
GZ cocycle on the periodic spectra; compare [32l Lemma 5.2], [T3 Proposition 3.2 
and Lemma 3.5]. The proofs would distract from the overall narrative flow of the 
paper; thus we postpone them until the end of the section. 
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Lemma 7.8. For each 6b > 7 r/ 4 , there exist constants m = m(6b), M = M[6b), 
and A = X{9b) such that the following holds. If db < da < 7r/2 and > X, then, for 
every z € (Tk, we have 

(7.18) < \x',{z)\ < 

where ^ = my/JI and S = My/Ji, as in the statement of Theorem [m 

Lemma 7.9. For all 9a and 9b as in the statement of Theorem \ 7.1\ there is a 
constant C > 0 such that for every n € Z with —2Fk < n < 2Fk with n ^ 0, one 
has 

(7.19) \\Z{n,0-,z)\\<C\n\^/^ 

for all z € ak, where t is as in (17.101) . and Z denotes the Gesztesy-Zinchenko 
cocycle associated to £q. Consequently, we have 

\\Z{n,m-,z)\\ < C^lnl^ 

for all z G (Tfe and all n,m with |n|, |m| < 2Fk and n 7 ^ 0 . 

Proof of Theorem [m Fix 9b and 9a £ {9b, tt/2) sufficiently large. The theorem 
consists of two nontrivial inequalities: a lower bound on /?“ and an upper bound 
on /?+. We will prove the lower bound using Theorem 12.11 and the upper bound 
using Theorem 12.51 Recall that Fk ^ 


Lower Bound. By Lemma lTTSl and the Mean Value Theorem, we have |i?fc | > 
for each band Bk of o-fc. Note that Fk ~ (p^. As a consequence of the 2Ffe-periodicity 
of the corresponding CMV operator and standard results in the theory of periodic 
operators (see, for example, [44l Theorem 11.1.1]), we know that ak consists of 2Fk 
nondegenerate closed intervals, known as bands. These facts imply that 

\<Jk\>8FkE-’^>F-\ 

where 77 = — 1, as in (17.101) . Note that the implicit constant on the right 

hand side depends on S, and hence on p and 9b. Combining this lower bound on 
the measure of the periodic spectra with the polynomial bounds from Lemma 17.91 
we get the desired conclusion from Theorem 12.11 by taking A/j = for each 

i? > 1, where i{R) G Z+ is the unique integer with < i? < F^. Specifically, 
with Rx = K tgtf ^ one has 


\Ar,\ = > R7 = 

which, by ( 12 .9L implies 

(l^l?„) (K) > \Ar^\K^ > 


Consequently, 




which proves the lower bound from (EH). 


Upper Bound. Define 


al := {z G C : \xk{z)\ < 1 + <5}. 
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Using the lower bound from ()7.18|) , one can argue as in [HI Proposition 3] and use 
the Koebe Distortion Theorem to see that for every d > 0 small enough, there is a 
constant Cs > 0 such that 


af CB(o,exp(c5r"/')) 

for all k G Z+. We now fix i5 > 0 small enough that has connected components 
for all n G Z_|_; the existence of such a S follows from the argument which proves 
[m Lemma 5]. Let > 0 be given, and put 


7 = 7(J^) 


log(^) 

2(1 + 1^) log((/3)' 


Thus, for all k G Z_|_ large enough (depending on i/), we have 


(7.20) ai U af+i C B (O, exp {CsF^)) . 

Let Kq be large enough that ze^o G at ^ whenever z G cr_i. Given K > Kq, 
define n = n{K) and N = N{K) by 


F'^ 1 F'>' 

<K< 


Cs 


Cs 


N = 2F, 


n-\-y/n'i 


where we abbreviate F^ := F^^^ for real numbers x. In particular, we have K ^ > 
CsF~'^, so, for complex numbers of the form with w G a-i, the escape 

condition (17.141) holds for some kg < n(K). On the other hand, (17.141) already 
holds with kg = S = 0 for all w G dlD \ cr_i by our choice of parameters, so, since 
escape is an open condition, we may enlarge Kg (if necessary) so that (j7.14l) still 
holds (with kg = S = 0) for all z of the form we^ with w G dll])\a-i and K > Kg. 
Consequently, for every 2 with \z\ = we have 

K{z)\>^l + -j 


for every integer q > n{K), by Lemma 17.61 For every e > 0, our choice of N{K) 
implies that there is a constant > 0 such that 

N{K) < 


for every K G Z_|_. Putting everything together, we have 


max 

0<n<CeA'T'"De 




> 


> 

xn(k) 





> 




for all 6 G [0, 27r), where we have used Lemma 17761 and N{K) — n{K) = ^yn{K) to 
obtain the final line. We then have 


max 

Q<n<c^K~i~^+‘ 




de 

27r 


< 



The right hand side decays faster than any negative power of A", and we can prove 
a similar bound on the other half-line using the reflection symmetry (17.11) . Thus, 
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we have < 7(u) ^ + e for all p > 0 by Theorem 12.51 Since this holds for all 

e > 0 and u > 0, we obtain /3^(p) < 7 ( 0 )“^, which is the upper bound in (17.91) . □ 

7.4. Band Combinatorics for Fibonacci Quantum Walks. Throughout this 
subsection, we adopt the standing assumptions Q < 9b < Oa < /2 and p > 32. 

Notice that the assumption /i > 32 implies ak-i Dcifc Hcrfe+i = 0 for all k > 0, since 
I{xk-i,Xk,Xk+i) > /i on (Tfc for all k. 

Let us say that a band Bk C cr^, is a type A band if Bk C and a type B 

band if B^ Q <Jk- 2 - The assumption on the 0’s means that ag consists of two bands 
of type A, and ai consists of four bands of type B. The first claim is an obvious 
consequence of 0 < 0b < 0a < t/2 , but the second requires an argument, which is 
supplied by the following lemma. 

Lemma 7.10. //O < 0b < 0a < 'x/2, then ct_i A ctq U cti. If, in addition, p > 32, 
then (To n (Ti = 0. 

Proof. As noted above, the assumptions on 6a and 0b immediately give ag C a-i. 
To prove the other inclusion, suppose z ^ cr_i, i.e., |Re(z)| > cos0b. Then the 
assumptions on the 0’s imply 

Re(z^) > cos(20b) > cos(0a + 0b) = cos 0a cos 0b — sin 0a sin 0b. 
Consequently, 

xi{z) = Re (z^) sec0a sec0b + tan0a tan0b 

> (cos 0a cos 0b — sin 9a sin 0b) sec 0a sec 9b + tan 9a tan 0b 

= 1 , 

so z ^ CTi. Thus (Ti C (T_i. With this in hand, we must have (Tq H cti = 0, since 
p>32. □ 

With the previous lemma in hand, we know that type A and B bands exhaust 
the periodic spectra at levels 0 and 1. To see that they exhaust the spectra at all 
levels, we use the following lemma, which is a CMV variant of [HU Lemma 5.3]. The 
proof is nearly identical to the proof in the Schrodinger case - we include detailed 
arguments for the reader’s convenience. 

Lemma 7.11. Assume 0 < 9b < 9a < Tr/2 and p > 32. For every k > 0: 

(1) Every type A band of au contains a type B band of ak +2 and no other bands 
of CTk+l or tTfc+2. 

(2) Every type B band of ak contains a type A band Bk+i C ak+i and two type 
B bands from of ak +2 whieh sandwieh Bk+i. 

Proof. Case A. Suppose Bk C ak is a type-A band. By definition, Bk C ak-i, so 
we must have Bk H Uk+i = 0, since p > 32; equivalently, jxfe+ij > 1 on Bk. There 
is a unique Zi € Bk such that Xk{zi) = 0. Using the trace map (17.41) . we have 

\xk+2{zi)\ = |a;fe_i(zi)| < 1, 

since zi € Bk Q Ok-i- In particular, BkC\ak +2 7 ^ 0- Notice also that when Xk = ±1, 
we have 

|a;fe+ 2 | > 2|a;fe+i| - \xk-i\ > 1, 

so any band of ak +2 which intersects Bk must be entirely contained within Bk. 
Finally, suppose Bk +2 is a band of ak +2 contained in Bk. By the IVT, Bk +2 
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contains a point Z2 for which Xk+2{z2) + Xk-i{z2) = 0 , so Xk{z2)xk+i{z2) = 0 by 
(d. Since Z2 ^ crfc+i, Xk{z2) = 0 . Thus, there is a unique band of crfc+2 contained 
in Bk- 


Case B. Suppose Bk C crj, is a type-B band, i.e. Bk Q o'k- 2 ] immediately, one has 
Bk n (Tk-i = 0. As above, if Xkizs) = 0, then 

|2;fe+l(^:3)| = \xk-2{z3)\ < 1, 

whence Bk fl ak+i 0- As in the previous case, \xk\ = 1 forces |a;fc+i| > 1, so any 
band of ak+i which meets Bk must be contained entirely within the interior of Bk- 
Moreover, running the argument above again, any band of crfc+i contained within 
Bk must contain a spectral parameter at which Xk vanishes, and so there is exactly 
one band of crfc+i contained in Bk- Now we consider bands of ak +2 which meet Bk- 
Notice that one may iterate the trace recursion to obtain 

(7.21) Xk +2 = 2xk+iXk - Xk-i = {4:xl - l)xk-i - 2xkXk-2- 

As before, any band of o'k +2 which meets Bk is completely contained in Bk ; specif¬ 
ically, if \xk\ = 1, then (17.2111 forces |a:fc+ 2 | >3 — 2 = 1. On the other hand, if 
Xk = ±1/2, then (17.211) forces |a;fe+ 2 | = |a;fc- 2 | < 1- Consequently, Bk contains at 
least two bands of ak +2 (note that we can’t have a single large band, becaue ak +2 
cannot meet ak+i, which contains the zero of Xk in Bk, which separates the points 
where Xk = ±1/2). It remains to see that we have no more than two bands of ak +2 
in Bk- Notice that 

{ 2 xk ± l)(a;fc+2 ± Xk-2) = (dXfc - l)(a;fe+i ± Xk-i). 

Now, fix a band Bk +2 ± Bk of < 7 ^+ 2 . We may choose a hxed sign so that a;fc+i±a;fc_i 
never vanishes on Bk+ 2 - By the IVT, there exists Z 4 € Bk +2 so that Xk+ 2 {zi) ± 
Xk- 2 {z 4 ) = 0. We must then have 1 ^( 24 ) = ±1/2. Since there are only two points 
in Bk where this happens, we are done - there cannot be more than two bands of 
ak+2 iri Bk- □ 

Lemma [7.111 implies that every band of ak is of type A or B for every fc > 0. 
Moreover, it also implies that ak consists of 2Fk disjoint bands for all k, that is, 
ak has no closed gaps. Finally, it also implies that (t_i A ak for all k, and hence 
(T_i D E by a strong approximation argument. 

7.5. Power-Lawr Bounds on the Growth of Fibonacci—Szego Cocycles. In 

this subsection, we prove Lemmas 17.81 and 17.91 Throughout, we have the standing 
assumptions that 7 r /4 < 9b < 9a < 7 r /2 and that ^ is large enough, where “large 
enough” is a 0 f,-dependent statement which becomes harder and harder to satisfy 
as 9b \ 7 r/ 4 . We do not attempt to find optimal constants, and instead opt for 
lucidity of presentation. Notice that 

+v^ + up' — 2uvw —1 = 1 

implies that 

w = uv± \/l -\- {1 — u^)(l — u^) 
via the quadratic formula. This motivates us to define 


g±{u, V, I) = uv ± + {1 — u2)(l — v^). 
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Lemma 7.12. For I > 4: and |m|, |t;| < 1, we have 


dg± 

du 


{u,v,I) 


dg± 

dv 


{u,v,I) 


< 1 . 


Moreover, 


^g^ / 


< 


1 

2V/' 


Proof. Since g±{u,v,I) = ^f±(2u,2v,2^/I) in Killip-Kiselev-Last’s notation, the 
first two inequalities are immediate from [321 Lemma 5.4] and the chain rule. In 
particular, note that / = A^/4 when passing between our notation and theirs. The 
third bound is obvious. □ 


Lemma 7.13. Suppose 'k 14 < Oh < tt/ 2 is given. There exist constants m = miOf), 
M = MiOb), and A = X{9b) such that if 6b < da < 7r/2 and fJ.> X, then the following 
statements hold with f = m^/Ji and S = M^/Ji. We have 

x'k+iiz) 

for all z € Bk+i if B^+i « type ^ 


for all z £ Bfe +2 if Bk +2 is a type B hand. 

Notation. We adopt the following notation, just over the course of the following 
proof. Given quantities / and g which depend on 9a and 9b (subject to all of our 
various standing assumptions thereupon), we denote / ~ (7 if for every 9b, there 
exists C = C{9b) with C~^f < g < Cf. 


hand, and 

4+2 (^) 


Proof. One can follow the proofs of [m Lemma 3.5] and |321 Lemma 5.5] and prove 
these statements by induction, but there are two new difficulties. First, the initial 
conditions are different in the CMV case, so we must verify that the induction is 
well-founded. Once this is done, we have to deal with extra terms which arise from 
the failure of I to be constant on the spectrum. To deal with the first base case, 
observe that 

x'oiz) _ sec 9a ^ 
x'_y(.z') sec 9b ^ 
for all z £ 9D. Next, for all z £ ai, we have 

— cos(0a — 9b) < Re( 2 ^) < cos(0a + 9b) < cos{29b) < 0 


by (17.121) and standard trigonometric identities. Consequently, |Im( 2 ;)| > sin^h and 
|Im( 2 :^)| > sin( 6 (Q — 9b), so we obtain 


x[{z) 


{2z — 2z ^)sec9asec9b 

x'-liz) 


(1 — z~^) sec9b 


Im(z^) 

Im(z) 

y/g- 


sec 9, 
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Now, we proceed inductively. Notice that the base case (and then each inductive 
step thereafter) implies \x').\ > 1 on Cfc as long as /r is large enough. Taking X[9b) 
sufficiently large, we have > 1 on cTfc for k = —1,0. 

In the argument below, we will need to estimate {dg±/dl)l'{z) for z G crk- Since 
\I'{z)\ = |Re(z)|K^ and |Im(z)| > sin0(, for all z S Ufe, we have 

< ^ M 

“ ‘iy/Jl sin^ 9b 2 sin^ 9b 

for all z G ak and all k > —1. 


(7.22) 


—I'(z) 
dl ^ ’ 


Case 1: Type A Band. Suppose Bfe+i is a type-A band of CTfc+i for some fc > 1; 
let Bk denote the band of ak with Bk+i C Bk- Consequently, we have BkHak-i = 0 
and Bk Q o'k- 2 , so Bk is a type B band. Differentiating the renormalization map 
(d and dividing by xj., we get 


''fc+i 


— ‘^^k—l 


‘2XkXf._^ Xi._2 


•^k -^k -^k 

To estimate the second term on the right hand side, note that Xk-i = 
g±{xk,xk- 2 ,l), so 


d9± , 


‘'k-l 


du 


‘' fc -2 


^9± , I dg± 
dv ^ dl ’ 


which implies 


‘'fc-i 


< 


‘' fc -2 


+ 1 + 




2|a;l.|sin^ 9b 


2/3 — 


<2 + 




2 sin^ 9b 


on Bfe+i, where we have used (I7.22L Bk+i ^ Bk ^ o-k- 2 , Lemma 17.121 and the 
inductive hypothesis. Consequently, 


(7.23) 


2XkXi^_^ Xf^_2 


<5 + 


sin^ 9b 


Since \xk\, |a;fe- 2 | < 1 on Bfc+i, one can use Xk-i = g±{xk,Xk- 2 ,I) to see that 

\/m - 1 < |a;fe-i| < v^+ 2 

on Bk+i- Combining this with ()7.23p and using the triangle inequality twice, we 
obtain 


V siv? 9b x^ V 


^ -x/m + 9 


sin 9b 


on Bk+i- Note that one needs 9b > 7r/4 for the lower bound to have any significance, 
even for large g. 


Case 2: Type B Band. Let us suppose that Bk +2 is a type-B band of <7^+2 for 
some fc > 0; Let Bk denote the band of ak which contains Bk+ 2 - There are two 
subcases to consider. (Notice that fc = 0 falls in Subcase 2ii, so one need not worry 
about the a;_ 2 ’s that might appear in Subcase 2i). 

Subcase 2i: Bk is type B. This implies Bk +2 ^ Bk C ak -2 and Bk fl ak-i = 0. 
As in the previous argument, we have 

y/g-l< \xk+i\ < v^ + 2 
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on Bk+ 2 , and we may estimate 
(7.24) 


^k-l 


< 2 




2 sin 9b 


on Bk+ 2 - Using the trace recursion (El and its derivatives, we obtain 
% = + 2xk-2 + - 1 )^ - 2xk^. 

^k 

This time, the Xk+i term dominates. Using (17.241) . the inductive hypothesis, and 
Bk +2 Q Bk C (Tfc _2 to estimate the remaining terms for z G Bk+ 2 , we have 


2xk-2 + {^xl-l)^-2xk^ 


<2+32+ 




2 sin^ 9b 


+ 2 


= 10 + 




2 sin^ 9b 


on Bk+ 2 - Using the triangle inequality again, we have 


4- 


2 sin 


^-14< 


Xi. 


fc +2 


< 4 


2 sin 9b 


y/Jk + 18 


on Bk+ 2 - 


Subcase 2ii: Bk is type A. We have Bk +2 ^ Bk ^ cTk-i, and Bk D ak -2 = 0- 
Differentiating the renormalization equation, we get 


^ = 2xk+i + 2xk 

^k 

Again, ^ - 1 < |a;fe+i| < + 2 on Bk+ 2 , 

Xk+i = g±{xk,xk-i,l), so 


<+i 4-1 



and, like the previous cases, we have 


Xk 


k+l 


on Bk+ 2 - Therefore, we have 


2xk ~ 


< 2 


+-1 




2 sin^ 9b 


<5 + 




sm 


on Bk+ 2 - Applying the triangle inequality two more times, we have 


2 - 


1 \ 


sin"^ 9b 


J 


y/JI-7 < 


'' fc +2 


< 2 


1 \ 


sin"^ 9b, 


+ 9 


on Bk+ 2 - 

Proof of Lemma \7^ This follows immediately from Lemma 1 7. 131 


□ 

□ 


Lastly, we will prove the power-law bounds on transfer matrices at spectral 
parameters in the periodic approximating spectra. To do this, we need trace bounds 
on these spectra, which the following lemma supplies. The key ingredient is that the 
escape condition is a bit simpler to check in our situation. Notice that something 
must be said in our situation, since some choices of 9a and 9b can lead to \x-i{z)\ > 1 
for z £ S, whereas one has a;_i = 1 in the Schrodinger case. 









































40 


DAVID DAMANIK, JAKE FILLMAN, AND DARREN C. ONG 


Lemma 7.14. If 0 < 9b < 9a < 1^12, > 32, and z € 9ID), the trace orbit 

(xk{z))k>-i is unbounded if and only if there exists fco > 0 such that 

(7.25) \xko{z)\, \xko+i{z)\ > 1. 

Consequently, if z € T, = Baa, one has |a;fc(z)| < k + 2 for all n > —1, where n is 
as in (17.811 . 


Proof. By Lemma |7.Ill we have cr-i 3 ak for all k, and cr_i 3 E. There are two 
possibilities: if z ^ cr_i, then one has z ^ T, and (I7.25|) holds for some fco. On 
the other hand, if z € cr_i, then |x_i(z)| < 1, and so the existence of kg which 
satisfies (17.251) is equivalent to the existence of kg which satisfies the general escape 
condition (17.141) . 

The claimed bound on |a:fc( 2 :)| for z € Baa follows from a simple observation. 
Namely, if \xk{z)\ > 1 for some z € Boa, then |a:fc±i(z)| < 1, whence 


where we have suppressed the dependence of the cc’s and I on z. 


□ 


Proof of Lemma \7.S\ Since this is nearly identical to m and m Lemma 3.2], we 
only sketch the main steps of the argument. The main idea is to begin by estimating 
the norms of the matrices M„, and then to decompose an arbitrary Z-type matrix 
over an even number of sites into a product of Mj’s by using the fact that every 
integer has a natural “base Fibonacci” expansion. Finally, one can interpolate to 
get transfer matrices over an odd number of sites. 

First, from Lemma 17.141 we have |xfe(z)| < k + 2 for every z G ak, where n is as 
in (Tra . As a first crude estimate, we have 

(7.26) ||M,(z)||<(k: + 2)^ 

for j < k and z £ ak- As in this follows from the recursive relationship among 
the M„-type matrices and the Cayley-Hamilton Theorem: 

~ Xj-lMj-2 — ^j-3- 

One also uses the fact that detM^ = 1 for all £, which implies HMjlgjl = ||Mj_ 3 ||. 
As in [23 Lemma 5] and [ITl Proposition 3.2], for every i > 2, j > 1, we may 


construct polynomials Pj^\ ■ ■ 

., Pj^^ of degree at most j in Xi-i,... ,Xj such that 


MkMk+j = PfM,+, + + Pfi. 

Moreover, 



(7.27) 

t 

{\xi-i\,... ,\x^+j\) < {2K + 5y, 


e=i 


where jPj denotes the polynomial obtained by replacing each coefficient of P by its 
absolute value. Indeed, (17.271) follows from a verbatim repition of the arguments 
of [13 [23 1 since all that is used in those references is Cayley-Hamilton and the 
trace bounds for indices bounded above by k. By repeating the next step of the 
argument, we obtain 


jjZ’(2m, 0; z)j| < for all 1 < m < Fk, 
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where d = {k + 2){2k + 5)^ and bm = 
Since Fn ^ we have 

bm 'F: 

for some constant D, whence 


max{n : F„ < m}; in particular, bm 


logm 

log(^ 


+ D 


< k. 


\\Z{2m,0-z)\\ < 


We obtain (17.191) for all n > 0 by interpolating to get odd values of n. The negative 
case follows from the reflection symmetry (EH. One may also obtain the same the 
result for T-type matrices by using (Iral) . □ 


8. The Thue-Morse Quantum Walk 

Let us consider a quantum walk as in Section [71 but for which the coins Ra and 
Rb are distributed according to an element of the Thue-Morse subshift, rather than 
the Fibonacci subshift. 

More precisely, consider the alphabet A = {a, b} as before. The Thue-Morse 
substitution is defined by 

S : a ab, b ^ ba. 

Iterating S on a, we obtain a sequence of words = S'"(a), and an infinite word 
mtm = lim Wn = abbabaabbaababba ..., 

n—^oo 

and we define the Thue-Morse subshift to be the space of all sequences over the 
alphabet A with the same local structure as utm, i.e., 

f^TM = {w G A^ : every finite subword of w occurs in utm} • 

As before, we can associate to each w G Utm a sequence of coins = Rt^r, and a 
sequence of Verblunsky coefficients as in Section [71 Now, consider the transfer 
matrices over the words More precisely, fix a = where Wq = MxmI^tm S 
rixMi and define 

M„( 2 ) = Z(2"+\0;z), n>0, . 2 GC. 

Also, define M_i as in the previous section. Let 

tn(z) = tr(M„(z)). 

The recursive relationship between the M„’s implies a recursive relationship 
amongst the t^’s, viz. 

(8.1) t„+i=t2_^(t„_2) + 2, n>2. 

The sets 

an = {z G dB : \tn{z) \ < 2} 

correspond to spectra of quantum walk update rules whose coins are distributed 
2"-periodicially according to Since the initial conditions are the same as in 
Section 13 we still have 

t-i{z) = Re(2;) secOb 
to{z) = Re(2) sec 0 a 

ti (z) = Re (z^) sec 6a sec 9b -I- tan 9a tan 9b 

Lemma 8.1. If zq G (?D satisfies tn(zo) = 2 for some n > 3 and t2{zo) 2 , then 
zq corresponds to a closed gap of an- 
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Proof. The proof is identical to the proof of the corresponding fact for 
Schrodinger/Jacobi transfer matrices, so we only sketch the details; compare [13 
Proposition 5.1]. Using the trace recurstion (EH), the stated conditions on zq im¬ 
ply that tfe(zo) = 0 for some 1 < k < n — 2. One can then use this to show that 
Mj^zq) = I for all j > fc -|- 1 (in particular, for j = n). Alternatively, one can 
differentiate 1° see that tj(zo) = 0 for all j > A: -|- 1. □ 

Theorem 8.2. For all ui G Otm; one has — 

Proof. This follows from Lemma 18.11 and Theorem 12.11 just as in the proof of 
Theorem 16.21 In particular, we take An = {zq} for all i? > 1, where zq G 9ID) 
satisifes t 2 {zo) ^ 2 and t„(zo) = 2 for some n > 3; compare [13 Theorem 4]. □ 

9. Ballistic Wavepacket Propagation for Periodic CMV Matrices 

We will say that a CMV matrix £ is p-periodic if = an for all n G Z. If 
necessary, double p so that it can be assumed to be even throughout this section, 
which is no problem, since nothing that we will do requires a minimal period. 
This can be thought of as a crystalline model, so the physical expectation is that 
resistance should go to zero and transport of wave packets should be free (ballistic). 
The following theorem establishes such a result 0 

Theorem 9.1. Suppose £ is p-periodic for some p G Z_|_. Then, for any ip G S, 
^^( 9 ) = ^^( 9 ) = 1 for all q>0. 

It is immediately clear that Theorem 19.11 generalizes to CMV matrices whose 
Verbhmsky coefficients are skew-periodic in the sense that an+p = wOn for some u! 
with |a;| = 1 and all n G Z. This also implies that a quantum walk with periodically 
distributed coins will exhibit ballistic transport. 

Corollary 9.2. Suppose £ is skew-periodic. Then, for any ip G S, = 

^^( 9 ) = 1 for all q>0. 

Corollary 9.3. Suppose U is the update rule of a quantum walk whose coins obey 
Qn-i-p = Qn for all nGZ. Then, for any ip G S, /3^(g) = P^{q) = 1 for all q > 0. 

Proof of Corollary 1 ,9. A Let w = fQj. some ^ G M, and introduce £ with coeffi¬ 
cients a defined by 

ctj = j G Z. 

It is not hard to see that a is honestly p-periodic, so we may apply Theorem l9.1l to 
deduce ballistic transport for £. Moreover, if we define the diagonal unitary T on 
e by 

7fc = exp (^^(-l)^+^ 

then £ = e~^^T*£V. Thus, we deduce ballistic transport for £. □ 

Proof of Corollaru \9.S\. Let A and £ = £a he as in Section [231 so that £ = A*UA. 
From the explicit form of A, we see that it maps S{£‘^{Z)) onto iS(£^(Z)(g)C^), and it 
does not affect transport exponents. Thus, it suffices to deduce ballistic transport 

8As we mentioned in the introduction, it is possible to derive this result from [T] Theorem 4]. 
We are grateful to Albert H. Werner for pointing this out to us. 


(p) , {5j,r6k) = Sj,kjk, 
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for £. Examining (I2.27|l . we see that an+ 2 p = oJotn for all n G Z, where w is given 
by 

sj= 

i=i 

Thus, S is skew-periodic, so ballistic transport for U follows from Corollarv l9.2l □ 


In order to apply the methods of [3] and m to deduce ballistic transport, we need 
to establish effective estimates on the Heisenberg evolution of the position operator, 
X, so we want to consider £X—X£ as in the other cases. We will actually work with 
a 2-block version of X, as this substantially simplifies the resulting commutators, 
without affecting the transport exponents. 

To set up notation, let = £^{Z) denote the space of finitely supported se¬ 
quences, and define a 2-block variant of the position operator X : by 


XSn = n5n := 2 


n 

L2J 


for n € Z. Clearly X is essentially self-adjoint on . For any observable (read: 
operator) A, denote its discrete Heisenberg evolution with respect to £ by 


A{j)=£^A£-\ jGZ. 

Obviously, is invariant under both £ and X, so that X(j) is a well-defined 
symmetric operator which is essentially self-adjoint on f' for all j G Z. As in m, 
the key ingredient is furnished by a variation on the theme of Asch-Knauf [3]. 


Theorem 9.4. There is a bounded, self-adjoint operator J such that 

s-lim —X(L) = J. 

L—^OO L 


Moreover, ker J = {0}^ and 


s-hm/(^ix(L)^ =/(J) 
for any bounded continuous function f. 

Remark 9.5. Notice that Theorem 19.41 does not follow from [151 Theo¬ 
rem 1.6(a),(b)] in any straightforward fashion. In order to deduce our theorem 
from theirs, one would need to be able to write £ = with H a Jacobi matrix (or 
some other self-adjoint operator of finite width), and it is not at all obvious that 
this is the case. Perhaps more importantly, even if such a proof could be found, it 
would add no insight to the unitary setting, whereas our proof explicitly constructs 
J out of computable functions of £. 

It is important to point out that the result is not trivial, even though we follow 
the outline of [15] . In particular, one needs an identity like (021) to make the 
argument work, and the more complicated structure of CMV matrices makes such 
an identity more difficult to prove. 
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9.1. Direct Integral Decomposition of Periodic CMV Matrices. The key 

ingredient here is to decompose a periodic CMV operator as a direct integral of 
unitary operators on C^. Since this is standard fare, we only summarize the the¬ 
orems and provide no proofs. See |441 Section 11.2] for more details on periodic 
CMV operators; see also [ 45 l Chapter 5] for the Jacobi version. 

Given a sequence a of Verblunsky coefficients, for notational convenience, intro¬ 
duce a, 6, c, and d as in [12], i.e., 

Un — 1 ; 17 — Pn^n—h dji — Tl € 

In terms of these parameters, the matrix representation of £ becomes 

/■. \ 


£ = 


ao 

^1 

di 



Co 

ai 

Cl 




h 

02 


ds 


C?2 

C2 

03 

C3 




64 

04 65 




d 4 

C4 as 


V 


■7 


If £ is p-periodic, then a mod-p variant of the Fourier transform diagonalizes £ 
into a multiplication operator on 'Hp = ([0, 27r),C^; ||), the space of functions 

/ : [0,27r) —>■ which are square integrable with respect to Lebesgue measure. 

More specifically, for each 0, denote 

\ 


£e = 


/ 

and define the mod-p Fourier transform T : P (Z) —> 'Hp by u i—>• u, where 


/ Oi 

^2 

Cl 

02 

^3 

ds 



dpC 

CpC *® 

d2 

C2 

O3 

C3 







64 

04 

bp -3 

1 

CO 





d4 

C4 

■ Op -3 

Cp —3 







■ bp-2 

Op -2 

bp—i 

dp—i 

Uie*® 

die*® 



dp-2 

Cp -2 

Qp—l 

bp 

Cp— 1 
dp 


.w-i: 




^-i£9 


1 < j < P, 6* e [0, 27r). 


Notice that our definition of £g has already implicitly used evenness of p; we will not 
comment on the parity of p again. The inverse of H is given by Hp B f f G 
where 


For any operator A on we denote by A its action on Fourier space, i.e., 

A = TAT~^. After simple calculations, we see that (^) = for each 

/ € Hp, and almost every 9 G [0,27r). It is typical (especially in the physics 
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literature) to view the Hilbert space T-Lp and the matrix £ as direct integrals of 
finite-dimensional objects: 


n 


/■® dO 

/ 



9.2. Proof of Ballistic Transport. The proofs of Theorem l9.1l and l9.4l are similar 
to those in m, so we will only comment in detail abont the challenges which 
one must overcome to rerun the Asch-Knauf-Damanik-Lukic-Yessen machine in 
the present setting, and we will be somewhat breezy with those details which are 
similar to the proofs in m- The main difficulty here is that time is discrete in our 
case, so we must use discrete derivatives with respect to time parameters, and this 
causes some minor headaches. 

Formally, let B = [£1, Y] = £X — X£. More precisely, B is initially defined on 
and then it is obviously a bounded operator thereupon (since £ is five-diagonal), 
and thus, B enjoys a unique extension to a bounded operator on t^(Z) by general 
nonsense. Now, observe that the Heisenberg evolution of B is (almost) the discrete 
derivative of the Heisenberg evolution of X: 


(9.1) 


X{k) = X + 



£*. 


The post-factor of £*, though unsightly, is ultimately irrelevant; more specifically, 
strong convergence of {B){L) clearly implies strong convergence of {B){L) ■ £* as 
i —>■ oo. A short calculation reveals 


/•. 


B = 


V 


\ 

0 0 2di 

0 0 2ci 

- 2&2 0 0 2^3 

-2^2 0 0 2c3 

-264 0 0 

-2di 0 0 


•7 


Since B is also clearly p-periodic, T also diagonalizes it into a direct integral as 
well. Specifically, we have (^) = Bgf{0) for all / S Bp and almost every 

0 G [0, 27r), where 


/ 

0 

2 ci 




-2dpe-*® 

0 \ 


-2&2 

0 

0 

2(^3 





2d2 

0 

0 

2 c3 








—2bp-2 

0 

0 

2dp—i 





—2dp-2 

0 

0 

2 cp_i 

V 

0 

2 die*® 




-2bp 

0 J 
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Now let Ve £ Cp^p denote the diagonal matrix with diagonal entries {ek^Vgtk) = 

= g2dfe/2jS/p^ J_g_^ 


/I 




o2iS/p 


o2iS/p 


J0{p-2)/p 


„i6l(p-2)/p 


„i0 


Define 


£0 = Vg-^SgVe, Bg = Vg-^BgVg 


By a direct calculation, one can check that 
(9.2) 

With this setup, we run the argument from m, which gives us 


p. = 

oO p 


L-l 


(9.3) 


lim — 
L —^00 L 


dXi 


^ £iBg£^^ = -ipY, 


t=Q 


J = 1 


for all 9 £ (0, 27r) \ {tt}, where Ai(0) < • • • < Am(0) are the distinct eigenvalues 
of £g and Pj{9) denotes projection onto the eigenspace corresponding to Xj{9). 
Notice that Pj{9) is a one-dimensional projection unless 0 = 0,7r. Let us briefly 
recapitulate the argument from |15j which proves (ESI). First, let Vj{0) be an 
analytic choice of a normalized element of ran(Pj(0)) = ker(fg — Xj{9)). As in |15) . 
with Vj{9) = V^^Vj{0), we have 

= ^{vj{0),£0Vji9)) 

= '-{v,{e),Bgv,ie)) 

= '-{v,{e),BgV,m- 

Additionally, 

A" I J2^0Bg£gWO)) = ^(y^ie),BgVkm 

_ f {vjiO),BgVki9)) a j = k 

1 0 otherwise 


Thus, (19.31) follows. 





SPREADING ESTIMATES FOR ID QUANTUM WALKS 


47 


Proof of Theorem \9.4\ Once one has (j9.3l) , then we can simply follow the proof of 
m Theorem 1.6]. In particular, we get the conclusions of Theorem 19.41 with 


J = —ipP 


' [0.277) 

More specifically, dominated convergence yields 




27T 


TZ\ 


(9.4) 


lim 

L—>-oo 


c ® ( 1 ^-1 

/ -E 

J [0,277) Q 


pin c-l\ 

£eBstg I — = -ip 


[0,277) 




i=i 


27r’ 


where the limit on the left hand side is taken in the operator norm topology. The 
left hand side is simply 

iL-l 


lim jE^W’ 

L—¥00 L/ 


^=0 


so, we can conjugate both sides of (19.41) by P to obtain 


L-l 


lim -y"B(i) = J£ 
L—^oo L ^ 


^=0 


The theorem then follows from (|9.1p and the arguments in [T^. 


□ 


Proof of Theorem \9.1\ This follows from Theorem 19.41 iust as in [15], since our 
choice of X does not change the transport exponents. □ 
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